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AUTOMATED PRELIMINARY DESIGN OF SIMPLIFIED WING STRUCTURES 


TO SATISFY STRENGTH AND FLUTTER REQUIREMENTS 

By W. Jefferson Stroud, Cornelia B. Dexter, 
and Manuel Stein 
Langley Research Center 

SUMMARY 

A simple structural model of an aircraft wing is used to show the effects of strength 
(stress) and flutter requirements on the design of minimum-weight aircraft-wing struc- 
tures. The wing is idealized as an isotropic sandwich plate with a variable cover thick- 
ness distribution and a variable depth between covers. Plate theory is used for the struc- 
tural analysis, and piston theory is used for the unsteady aerodynamics in the flutter 
analysis. Mathematical programming techniques are used to find the minimum-weight 
cover thickness distribution which satisfies flutter, strength, and minimum-gage con- 
straints. This report presents the method of solution, some sample results, and the com- 
puter program used to obtain these results. 

The results indicate that the cover thickness distribution obtained when designing 
for the strength requirement alone may be quite different from the cover thickness distri- 
bution obtained when designing for either the flutter requirement alone or for both the 
strength and flutter requirements concurrently. This conclusion emphasizes the need for 
designing for both flutter and strength from the outset. 

INTRODUCTION 

The speed and storage capacity of modern-day computers, together with the devel- 
opment of efficient computer-oriented decision-making strategies, have encouraged 
researchers to seek ways to use the computer for not only analysis purposes but also for 
portions of the design process. Research is underway in this area because automated 
design methods have the potential for providing superior designs in a shorter time and at 
a lower cost than is possible with more conventional design techniques. One of the objec- 
tives of automated design is to incorporate interdisciplinary design requirements into the 
design process. The present report considers one example of an interdisciplinary design 
problem - that of designing a wing structure to meet a strength (stress level) require- 
ment and a flutter requirement. In this example, an optimization procedure is used to 
determine the lightest structure which will meet these requirements. 



The purpose of this report is to show the effects of strength (stress) and flutter 
requirements on the design of simplified aircraft- wing structures. The wing is idealized 
as an isotropic sandwich plate with a variable cover thickness distribution and a variable 
depth between covers. Plate theory is used for the structural analysis, and piston theory 
is used for the unsteady aerodynamics in the flutter analysis. Mathematical programming 
techniques are used to find the minimum -weight cover thickness distribution which satis- 
fies flutter, strength, and minimum-gage constraints. This report presents the method 
of solution, some sample results, and the computer program used to obtain these results. 

Some previous work dealing with the design of aircraft structures to satisfy aero- 
elastic constraints may be found in references 1 to 5. Some studies of automated struc- 
tural design for strength, buckling, dynamic and/or deflection constraints are contained 
in references 6 to 12. 


SYMBOLS 

The physical quantities used in this paper are given in both the International System 
of Units (SI) (ref. 13) and in the U.S. Customary System of Units. The measurements and 
calculations were made in U.S. Customary Units. Appendix A presents factors relating 
these two systems of units. 

a speed of sound 

C design variables 

c local chord length 

D local flexural stiffness of wing 

E modulus of elasticity 

h total local depth of wing 

M Mach number 

m mass of wing per unit area 

N number of spanwise stations on the wing (see fig. 3) 

P sum of the objective function and the penalty function (see eq. (16)) 

p normal loading per unit area 
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r penalty -function weighting factor (see eq. (16)) 

t local thickness of wing covers 

v velocity 

v £ flutter velocity 

W, w assumed normal deflection of wing neutral surface (see eqs. (1) and (2)) 

X, Y,Z chordwise, spanwise, and normal axes, respectively 

x,y chordwise and spanwise coordinates, respectively 

x^ value of x at wing leading edge 

x,p value of x at wing trailing edge (zero for all calculations presented herein) 

y g value of y at wing tip, semispan 

a real part of flutter determinant 

/3 imaginary part of flutter determinant 

y ratio of specific heats 

e spanwise finite -difference increment 

e x ,6y.r ^ components of normal and shearing strain in wing cover plates 
j u Poisson's ratio 

p mass density of air 

o cr critical value of the stress defined by the Von Mises stress condition 

°x>°y’ T xy components of normal and shearing stress in wing cover plates 
t time 
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<p^ deflection coefficients (see eq. (2)) 

co frequency, rad/s 

oo ^ flutter frequency 

Dots indicate differentiation with respect to time. 

ANALYSIS 


The analysis procedures used in this report are plate theory for the wing structural 
analysis, piston theory for the unsteady aerodynamics in the flutter analysis, and mathe- 
matical programming techniques for finding the minimum-weight design. These three 
procedures are discussed briefly in the next three sections. Additional discussion and 
detailed definitions of matrix expressions are given in appendix B. 


Plate Theory for Wing Structural Analysis 

The wing is assumed to be sufficiently thin to be analyzed with plate theory. It is 
further assumed that internal bulkheads, ribs, and spars exist which prevent transverse 
shear deflections, but any other load- carrying capacity of these structural elements is 
neglected. All load is assumed to be carried by a set of wing cover plates — an upper 
cover plate and a lower cover plate. The wing structural analysis is concerned only with 
these cover plates. The wing is assumed to have a symmetric biconvex airfoil shape and, 
therefore, the distance between covers varies over the surface of the wing (fig. 1). The 
covers have a variable thickness distribution which is adjusted in the design procedure to 
provide the minimum-weight design. 



Section A A 

Figure 1.- Wing coordinate system. 
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The wing structural analysis is based upon the approach presented in reference 14. 
The normal deflection W of the wing neutral surface is assumed to be given by 

W(x,y,r) = w(x,y)e ia,T (1) 

where 

w(x,y) = <p Q (y) + + x 2 cp 2 ( y) (2) 

The assumed solution includes terms only to the second degree in the x (chordwise) 
direction. Finite differences in the y (spanwise) direction are used to solve for the 
unknowns <Pq (y), <pj(y), and <p 2 (y). The wing coordinate system and deflection shape 

are shown in figures 1 and 2, respectively. The location of the Y-axis is arbitrary. 



Figure 2.- Diagram of undeflected and deflected neutral surface of 
wing showing assumed parabolic chordwise-deflection shape. 


An energy approach (ref. 14) is used to obtain the equilibrium equations for a 
variable-thickness isotropic plate or for an isotropic built-up plate with variable depth. 
The matrix form of these equilibrium equations is 


A 00 A 01 A 02 


<Pq 


PO 

A 10 A 11 A 12 


*1 

= € 4 

Pi 

& 

o 

>— * 

& 

i 


cp 2 


p 2 


( 3 ) 
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in which the matrix is composed of terms involving the stiffness and derivatives 

of the stiffness, the vector describes the deflections, the vector j’p.j is the nor- 

mal loading, and e is the spanwise finite difference increment shown in figure 3. The 
matrix expressions are presented in detail in appendix B. The normal loading used in 
the strength analysis is a uniform normal loading. The wing is assumed to be fixed at 
the root. 


Station 

number 



Figure 3<- Location of stations on planform of half -wing. 

Once the deflections have been calculated from equations (3) and (2), the strains in 
the wing covers are calculated by using the following strain-displacement relations: 


£ h 9 2 W 

(4) 

h a 2 w 

V ^ 2 3y2 

(5) 

. 9 2 W 

y = q=h 

xy 9x 9y 

(6) 


where the upper and lower signs are for the upper and lower wing covers, respectively, 
and h is the local distance between the upper wing- cover midplane and the lower wing- 
cover midplane. The stresses can then be calculated by using the following stress-strain 
relations: 



(7) 
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(8) 


°y = * M£x ) 




xy 2(1 + m) *y 


(9) 


These stresses are used in the strength constraint discussed in a subsequent section 
entitled "Mathematical Programming Techniques for Determining Minimum-Weight 
Design." 


Piston Theory for Unsteady Aerodynamics in Flutter Analysis 

In contrast to the uniform normal loading used in the strength analysis of the wing, 
the normal loading used in the flutter analysis is a function of x, y, and t (time). 

The total loading that acts on an airfoil in motion is the sum of the inertia and unsteady 
aerodynamic loadings. (Structural damping is neglected.) Based on second-order piston 
theory for symmetrical airfoils (refs. 15 and 16), the sum of these two loadings is 

p(x,y,T) = -mW - 2pa^l + ^+1 M + v -^jw (10) 

The unsteady loading p(x,y ,r) can be taken to the left side of the equilibrium equation 
(eq. (3)), thereby allowing the equation to be written in the form 


o 

o 

> 

o 

>-* 

A 02 


B oo 

B 01 

> 

B 02 


% 

> 

o 

A 11 

A 12 

+ e4 

B 10 

B 11 

B 12 

> 

*1 

A 20 

A 21 

^2 


B 20 

B 2 1 

B 22 


cp 2 
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-> 
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Detailed expressions which define the components of equation (11) are given in appendix B. 

The flutter speed Vf and the flutter frequency are calculated by setting the 
complex determinant of the matrix of coefficients of <p ^ in equation (11) equal to zero — 
that is, 



a + i0 = 0 


( 12 ) 
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There are, therefore, two equations 


0=0 

(13) 

0 = 0 

(14) 


in the two unknowns Vf and First, the solution neighborhood for equation (12) is 

established by evaluating the determinant for many combinations of co and v. The 
Newton-Raphson method is then used to determine u>f and Vf precisely. An example 
of this solution technique is given in appendix C. 

Mathematical Programming Techniques for 
Determining Minimum-Weight Design 

The thickness distribution of the wing covers is adjusted to give the minimum- 
weight design which satisfies all the design requirements. The thickness distribution is 
assumed to be of the form 

t(x,y) = ^ Y CijxV (15) 

i j 


The coefficients Cjj are design variables which are adjusted to give the minimum- 

weight design. Henceforth, the design variables will be referred to as the vector of 
design variables denoted by C. The detailed expressions used for t(x,y) are given in 
a subsequent section entitled "Applications and Results.” 

In this report the constraints, which are the design requirements, are introduced 
into the problem by way of an interior penalty function. This method of treating con- 
straints is sometimes called "SUMT" — Sequential Unconstrained Minimization Technique 
(refs. 17 and 18). In SUMT, the objective function, which is taken to be the weight in this 
study, is augmented by a penalty term that accounts for the constraints. The combination 
of these two terms is called the P function: 



= Weight! C ) +■ r 


i S A C 


(16) 


The term 



is the penalty term in which 


r 


is a weighting factor and 
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are the inequality constraints. The inequality constraints are formulated in such a way 
that they may be written as 

gj(c) > 0 (17) 

The penalty term repels the P function from the inequality constraints and trans- 
forms the problem from a constrained minimization problem to an unconstrained problem. 
By carrying out a sequence of minimizations of the P function corresponding to a 
sequence of reductions in the value of the weighting factor r, it is possible to obtain the 
vector of design variables which satisfies the constraints and gives the minimum value of 
the objective function. The constraints used herein are as follows: 

(1) For flutter, the flutter speed Vf shall not be less than some prescribed critical 
flutter speed Vj cr . This constraint can be expressed mathematically as 

< 18 > 


(2) For strength, at no point on the wing covers shall the stress exceed the pre- 
scribed critical stress, or 


^2 


A J J „ cr 


> 0 


"cr 


where A is the wing area, 


CT= K 2 -°x a y + 0 y 2 + 3 T xy 2 


1/2 


(19) 


which is the Von Mises stress condition, and cr cr is the critical value of cr. 

(3) For minimum gage, at no point on the wing covers shall the thickness t be less 
than some prescribed minimum thickness t m i n , or 


^3 


iff 


> 0 


dA 


A 1 


t 


min 


( 20 ) 
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Once a set of design variables has been defined (eq. (15)), and the design require- 
ments have been formulated as constraints (eqs. (18) to (20)), the wing can be designed 
by using mathematical programming methods. Additional comments on the math pro- 
gramming approach used are given in appendix B. 

A computer program denoted "SWIFT" which incorporates the structural analysis, 
flutter analysis, and optimization techniques just discussed is presented in appendix C. 

In order to check the validity of the structural analysis, natural -vibration frequen- 
cies for a trapezoidal flat plate were calculated by using the preceding analysis and are 
compared with experimentally determined natural frequencies in appendix D. A theoret- 
ical flutter speed and frequency for that plate are also presented in appendix D but are 
not compared with experimental results. 

APPLICATIONS AND RESULTS 

In order to show the effects of flutter and strength requirements on the design of 
aircraft wings, these requirements were used in the calculation of several minimum- 
weight wing cover thickness distributions. In these sample calculations, the number 
and choice of design variables is varied and the design requirements are varied. The 
wing properties used in all calculations are presented in figure 4. For the examples 
presented herein, the Y-axis is chosen to lie along the wing trailing edge. The number 
of spanwise stations N used in the structural analysis is six. 

The following specific examples were considered: 

First, the following set of design requirements was chosen: 

(1) For flutter, the wing is to be flutter-free at a dynamic pressure of 
160 kN/m 2 (3340 lbf/ft 2 ) at an altitude of 7620 meters (25 000 ft). (This cor- 
responds to M = 2.46.) 

(2) For strength, the wing is to support a uniformly distributed loading of 
6.89 kN/m2 (144 lbf/ft 2 ). 

(3) The minimum-gage requirement is 0.0508 centimeter (0.0200 in.). 

These design requirements, together with a sixteen-design-variable polynomial, were 
used to obtain a flutter design (requirements 1 and 3), a strength design (requirements 2 
and 3), and a design (denoted as a "strength-flutter design" or "combined design") which 
meets the flutter, strength, and minimum-gage requirements. 

Second, these same design requirements, together with a six-design-variable poly- 
nomial, were also used to obtain a flutter design, a strength design, and a strength-flutter 
design. 
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X 



Figure 4.- Wing description and cover -material properties used 
in the sample calculations. Fuel mass, I 50 000 lbm (68.0 Mg). 
Properties of titanium material: E = 16 400 ksi (113 GN/m 2 ); 

p = 0.3; o cr = 125 ksi (863 MN/m 2 ) • density, 0.160 lbm/in3 

(4430 kg/m3). 


Third, the flutter and strength design requirements were varied and the mass of 
minimum -weight wings designed by using the six-design-variable polynomial was calcu- 
lated as a function of the design requirements. 


Results Obtained With Sixteen Design Variables 

The thickness distribution chosen for the first example is a sixteen-term polyno- 
mial. That polynomial is 


t = Cl + C 2 |l-fVc 3 (l-ff + c 4 (l-l -' 3 


° s +c « f 1 '^) + ° 8 


c 9 + C lo( 1 " + C llf 1 " +C 12 [ 1_ “ 


1--^ 
V X L 


ya 


C13 + C - t —\ + C-i q(l - — 1 + Ci ft [1 - — 


'15 


ys 


'16 


ys 


1 - 
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2x 


2x 
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in which x is the chordwise coordinate, y is the spanwise coordinate, is the 
value of x at the wing leading edge, yg is the semispan, and Cj (i = 1,2,3,. . .,16) 
are the design variables. For clarity, the spanwise and chordwise variations corre- 
sponding to each term in equation (21) are shown in figure 5. The thickness distribution 
for the upper cover is the same as for the lower cover. 



“13 1 


+ C, 


+ C 


15 L 


+ C 


16 L 




— tfd 1 


Figure 5 • - Spanwise and chordwise variations for each term in equation (21). 
R = root; T = tip; = leading edge; = trailing edge. 


Contour plots showing wing cover thickness distributions for the flutter, strength, 
and combined designs are shown in figure 6. These plots show contours for a single 
cover. If only the flutter and minimum-gage requirements are considered, the lightest 
cover thickness distribution that can be obtained from the sixteen design variables 
described in equation (21) is shown in figure 6(a). The mass of the semispan wing (that 
is, the mass of both the top and the bottom cover) is 3160 kilograms (6960 lbm). This 
flutter design will support a uniform static pressure loading of only 2.82 kN/m^ 

(58.9 lbf/ft2). If only the strength and minimum-gage requirements are considered, the 
lightest cover thickness distribution for the sixteen prescribed design variables is shown 
in figure 6(b). The cover mass is 3220 kilograms (7100 lbm). This strength design will 
flutter at a dynamic pressure of 51.2 kN/m2 (1070 lbf/ft^). The cover thickness distri- 
bution for the flutter design is quite different from that for the strength design. The 
design which satisfies both the flutter and strength requirement is shown in figure 6(c). 
The cover mass for this design is 4030 kilograms (8890 lbm). 


Because of the large mass of fuel (68.0 megagrams (150 000 lbm)) carried in the 
wing, the contribution of the wing covers to the total mass distribution is small but is, 
nevertheless, included in the flutter analysis used to obtain these results. The primary 
effect of the thickness distribution is to establish the stiffness distribution. 


12 



(. 076 ) (. 076 ) ;(. 15 ) 



(. 15 ) 1 . 076 ) (. 076 ) 1 . 15 ) (. 23 ) (. 23 ) 




Figure 6.- Contour plots showing wing cover thickness 
distributions for flutter, strength, and combined 
designs obtained by using sixteen design variables. 
Contour labels are in inches (centimeters). 
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Results Obtained With Six Design Variables 

Calculations were also made with a judicious combination of certain terms in equa- 
tion (21). The thickness expression which was used is 



in which t m ^ n = 0.0508 centimeter (0.0200 in.), the minimum-gage requirement. The 
spanwise and chordwise variations for each term in equation (22) are shown in figure 7. 
The thickness distribution described by equation (22) is minimum gage along the leading 
and trailing edges and at the tip - a reasonable limitation. 


R T R T R T 1 f X L X T 



Figure 7 .- Spanwise and chordwise variations for each term in equation (22). 
R = root; T = tip; x L = leading edge; x T = trailing edge. 


All of the thickness variations described by equation (22) can also be obtained from 
equation (21). Therefore, the design obtained by using equation (22) cannot be lighter than 
the design obtained by using equation (21). Equation (22), however, contains only six 
design variables, compared with sixteen design variables in equation (21). Since com- 
puter time for the design process depends strongly upon the number of design variables, 
it is advantageous to use a small number of well-chosen design variables rather than a 
large number of arbitrarily chosen design variables. 

The design requirements used for this second example are the same as those pre- 
viously stated. Contour plots showing wing cover thickness distributions for the flutter, 
strength, and combined designs are shown in figures 8(a), 8(b), and 8(c), respectively. As 
in the preceding example, the cover thickness distribution for the flutter design is quite 
different from that for the strength design. 
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Figure 8.- Contour plots showing wing cover thickness 
distributions for flutter, strength, and combined 
designs obtained by using six design variables. 

Contour labels are in inches (centimeters). 

For the flutter design, the new six-design-variable thickness distribution is 10 per- 
cent heavier than the sixteen- design- variable thickness distribution. The new strength 
design is only 5 percent heavier than the sixteen-design-variable case. The new combined 
design is only 2 percent heavier than the sixteen-design-variable case. 

Results Obtained for Various Strength and Flutter Design Requirements 

The wing properties given in figure 4 were used in calculating several other 
minimum -weight thickness distributions for various strength and flutter design require- 
ments. All of these calculations were made with the six-design-variable polynomial 
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given in equation (22). The results are presented in figure 9 in which minimum -weight 
wing cover mass is shown as a function of flutter (abscissa) and strength (ordinate) 
design requirements. Lines of constant cover mass are shown as dashed lines. The 
solid curve which is almost horizontal represents designs dictated by flutter require- 
ments, and the solid curve that is almost vertical represents designs dictated by strength 
requirements. Points a, b, and c indicate the flutter, strength, and combined designs 
shown in figures 8(a), 8(b), and 8(c), respectively. 


Design flutter dynamic pressure, 


kN/m 2 


0 


100 200 300 



Design static wing 
2 

loading, kN/m 


Figure 9.- Minimum-weight wing cover mass as a function of design flutter 
dynamic pressure (flutter requirement) and design static wing loading 
(strength requirement). 
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The results presented in figure 9 indicate that for this example there is little inter- 
action between strength and flutter design requirements. For a strength-designed wing, 
improving the strength in the best possible way improves the flutter characteristics only 
slightly. For a flutter-designed wing, increasing the flutter speed in the best possible 
way has a negligible effect on the strength. The fact that designing for strength is quite 
different from designing for flutter emphasizes the need for designing for both strength 
and flutter from the outset. 

There are many ways to obtain a design which satisfies both the flutter and strength 
design requirements given by point c in figure 9. One approach is to start with the 
strength design shown at point b and increase the thickness uniformly over the entire wing 
until the flutter requirement is satisfied. For the present set of design conditions, the 
resulting design is 38 percent heavier than the minimum-weight combined design at 
point c. 

A second approach is to use the envelope of the thickness distributions for the flut- 
ter and strength cases. That is, at each point on the wing the cover thickness is taken to 
be the larger of the two following thicknesses: (1) the thickness at that point for the flut- 
ter design, and (2) the thickness at that point for the strength design. In the present case, 
the resulting envelope design does meet the design requirements and is only about 3 per- 
cent heavier than the combined design at point c. However, all envelope designs do not 
meet the design requirements. For example, if the strength design requirement given by 
point c in figure 9 is reduced 50 percent to 3.45 kN/m 2 (72 lbf/ft 2 ) without changing the 
flutter design requirement, the envelope design corresponding to the new set of design 
requirements is about 2 percent lighter than the minimum -weight combined design for that 
same set of design requirements. However, the envelope design flutters at a dynamic 
pressure about 9 percent lower than the flutter design requirement. 

A third possible approach is to use the strength design as a minimum-gage require- 
ment and, by using an optimization technique, add material to meet the flutter require- 
ment. Such a design would be lighter than an acceptable envelope design. There are, of 
course, many other approaches. 

The intent of this report is not to prove that a combined design approach as char- 
acterized by the math programming techniques used herein is superior to any other design 
approach. There are, however, two items to consider with respect to the combined design 
approach. First, for the set of design variables considered, the combined design approach 
will give the lightest design that meets all design requirements. In this simple model 
which combines two disciplines - structural analysis and aeroelasticity — the weight 
penalty caused by designing separately and sequentially rather than concurrently may be 
small. However, as the model becomes more complex and other disciplines are added, 
that penalty may increase. Second, for the optimization techniques described herein it is 
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faster to carry out a combined design than to design for strength and then to design for 
flutter. It is realized, however, that there are specialized design techniques that are 
tailored for a particular discipline. These techniques may be well-suited for sequential 
design and may provide an adequate design with shorter computer run times. The effec- 
tiveness of these design techniques can be assessed by comparing the results obtained by 
using these techniques with the results obtained by using a combined-design math- 
programming approach. 


CONCLUDING REMARKS 

A simple structural model of an aircraft wing has been used to show the effects of 
strength and flutter requirements on the design of minimum -weight wing structures. The 
wing structural analysis is based on plate theory. Piston theory is used in the flutter 
analysis. This report has presented the method of solution, some sample. results, and the 
computer program used to obtain these results. 

Wings were designed by adjusting the coefficients of polynomials which describe the 
wing cover thickness distribution. Sample calculations were carried out for a polynomial 
containing six coefficients (design variables) and for a polynomial containing sixteen coef- 
ficients. Mathematical programming techniques were used to determine the values of the 
coefficients which gave the lightest cover thickness distribution that satisfied flutter, 
strength, and minimum -gage requirements. 

The results indicate that the cover thickness distribution obtained when designing 
for the strength requirement alone may be quite different from the cover thickness dis- 
tribution obtained when designing for either the flutter requirement alone or for both the 
strength and flutter requirements concurrently. This conclusion emphasizes the need 
for designing for both flutter and strength from the outset. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 27, 1971. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General Confer- 
ence on Weights and Measures, Paris, October 1960 (ref. 13). Conversion factors for the 
units used herein are given in the following table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 
(**) 

Force 

lbf 

4.44822 

newtons (N) 

Length 

f ft 

0.3048 

0.0254 

^ meters (m) 

Mass 

lbm 

0.45359 

kilograms (kg) 

Stress, modulus . . . 

ksi 

6.895 x 10 6 

newtons/meter2 (N/m2) 

Density 

lbm/in3 

27.68 x 10 6 

grams/meter 3 (g/m3) 

Pressure 

lbf/ft2 

47.88 

newtons/meter2 (N/m2) 

Speed 

ft/s 

0.3048 

meters/second (m/s) 


^Multiply value given in U.S. Customary Units by conversion factor to obtain 


equivalent value in SI Units. 

** Prefixes to indicate multiples of units are as follows: 


Prefix 

Multiple 

giga (G) 

10 9 

mega (M) 

10 6 

kilo (k) 

10 3 

deci (d) 

io- 1 

centi (c) 

lO- 2 

milli (m) 

CO 

1 

o 

t-H 
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APPENDIX B 


DEFINITIONS OF MATRIX EXPRESSIONS AND FURTHER 
DISCUSSION OF ANALYSIS 

Plate Theory for Wing Structural Analysis 
For convenience, the plate equilibrium equation (eq. (3)) is repeated here 


A C 

A 1 

Aq 


A e 

A 1 

Ar 


A C 

A 1 

Ac 


2 

<?0 


Po 

2 

n 

= € 4 

Pi 

2 

<P2 


P2 


The matrix components of equation (Bl) are 


a 0 

D 0 
~ — 

— — 
a l 

1 O 1 
1 0 , 


[ago] = [do 

H - [°0 

[a 02 ] = j^o ^ a 2 

[An] = Jdo 

N - [°0 

[a 22 ] = [do 




D f 


+ 2 1 ±e* 


ao 


H 

[do] 

+ 2(1 - p)e 2 




S 

[W 

T 

+ 4(1 - n)e 2 

fa] 



1 Tjl ) 
1 & 1 


f 

+ 8(1 - p)e 2 

D 3 

H 

[”•] 


+ 2/ne 2 
+ 2pe 2 


D 


a l 

~\ 

D 0 


a 2 


D 0 a 2 ^ j + 4e 4 a Q 


[aio] = [ao] [a2o| = [a-os] [A21] = [ai 2 ] 


J 


(Bl) 


(B2) 


Primes indicate the transpose of a matrix, and e is the finite-difference increment in 
the spanwise direction. The finite-difference increments are shown in figure 3. 
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In this report the wing is assumed to be fixed at the root. Because more general 
boundary conditions are considered in reference 14, some of the expressions used herein 
differ slightly from the comparable expressions in reference 14. For the most part, 
however, the quantities defined herein are the same as those defined in reference 14. 

The matrices in equations (B2) are defined as follows (where the number of rows and 
columns, in that order, are indicated below the equation numbers): 


2 


-2 1 
1 -2 



1 


1 

-2 1 

• • • 

1 -2 1 

1 -2 
1 


(B3) 

(N+1)X(N+1) 


1 


-1 

1 



-1 

1 


-1 



(B4) 

(N+l)xN 



k=0,l,2,3,4 




a k,N-l 


2 a k,N 


(B5) 

(N+1)X(N+1) 
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in which 


0 


L a kJ 

k=0,l,2 


a k,l 

a k,2 


a k,N-l 

| a k,N 0 


(B6) 

(N+l)x(N+l) 



(B7) 

NXN 


(B8) 

(N+l)x(N+l) 



Eth 2 ) 2Et 3 
2(l - m 2 ) 12(l - ix 2 ) 


(B9) 


(BIO) 
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and 

E modulus of elasticity 

t local thickness of either the upper or lower cover (airfoils are symmetrical) 

h total local depth of wing 

H Poisson's ratio 

The quantity a^ g is the value of a k at the root, y = 0; a, x is a^ at y =■! . e; 

2 

a^ j is a^ at y = 1 • e; and a^ ^ is a^ at the tip. The wing tip is located at 
y = N • e. (See fig. 3.) Several of the quantities which follow are defined in the same way. 

The column vectors in equation (Bl) describe the displacements and loading and are 
given by 



and 



(Bll) 


(B12) 


23 



APPENDIX B — Continued 


in which 



(B13) 


This definition of p k differs from that given in reference 11. As a result, there is a 

factor in equation (Bl) compared with a factor in the equivalent equation 
(eq. (30)) of reference 14. The quantity p in equation (B13) is the intensity of the nor- 
mal loading and has units of force per unit area. 

In the present report the static normal loading used in the strength analysis of the 
wing is assumed to be a uniform loading; therefore, 


p k (y) = p J 


<r(y) 


dx 


x T.(y) 


(B14) 


This completes the definitions of terms that appear in equation (Bl). The deflec- 
tions are calculated from equations (Bl) and (2), and the strains are calculated from 
derivatives of the deflections (eqs. (4) to (6)). The stresses are calculated from the 
stress-strain relations given in equations (7) to (9). 

Piston Theory for Unsteady Aerodynamics in Flutter Analysis 

When the expression for W(x,y,r) given in equation (1) is substituted into equa- 
tion (10) the following equation results: 


P(x,y,T) 


mo)2w - 2pa 1 +^~ 

- M— ^ 

(■ 8w\ 

iww + v — 

\ 4 

dx) 

\ ax/ 


e iwt 


where the positive direction is upward and 
m mass of wing per unit area 

p mass density of air 

a speed of sound 

y ratio of specific heats 


(B15) 
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M Mach number, v/a 

h total local depth of wing 


v velocity 


As noted in the body of the report (see eq. (11)), the unsteady loading p(x,y,T) can be 
taken to the left side of the equilibrium equation (eq. (3) or (Bl)), thereby allowing the 
equation to be written in the form 


AqO ^1 ^2 B 00 B 01 B 02 h> 


a 20 A 21 B 20 b 21 b 22 Wl 


in which the components for Ajj are given in equation (B2) and the components for Bjj 
are as follows: 
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’k,2 


LM 

k=0,l,2,3,4 


k,N-1 


1 s 

2 b k,N 


(B18) 

(N+l)x(N+l) 


k,l 


k,2 


1 


k=0,l,2,3,4 


k,N-l 


1 T? 

2 k,N 


(B19) 

(N+l)x(N+l) 


H, 


M 


H 


k,2 


L Hk J 

k=0,l,2,3,4 


H 


k,N~l 


1 H 

2 k,N 


(B20) 

(N+l)x(N+l) 



(B21) 


(B22) 


(B23) 
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This completes the definition of terms that appear in equation (B16). The determi- 
nant of the matrix of coefficients of <pj in equation (B16) is the flutter determinant 
(eq. ( 12 )). 


Mathematical Programming Techniques for Determining 


Minimum -Weight Design 


In the interior penalty-function method denoted as "SUMT," the vector of design 
variables corresponding to the minimum value of the objective function is found by 
carrying out a sequence of minimizations of the P function. The P function P^C,r k 
is defined as 



(B24) 


In equation (B24) , 



is the objective function, C 


is the vector of design variables, 


and 



is the penalty term in which 



are the inequality constraints 


and r k is the constraint weighting factor. The subscript k placed on r refers to 
the sequence of values that r assumes during the minimization process. Let Cj be 
the vector of design variables which satisfies the inequality constraints and minimizes 
the P function for an initial value of r. Let that value of r be denoted as r j. When 
Cj is found, the value of r is reduced to ^2 < r^ and a search is carried out to find 
C 2 , which is the vector of design variables that minimizes pfc,^). The starting point 
for that search is Cj. When C 2 is found, the value of r is reduced to rg < r 2 < r j. 
Starting with C 2 , a search is carried out for C 3 . This sequence of searches is con- 


tinued until the penalty term contributes little to the P function which is then approxi- 
mately equal to the objective function. 


The manner in which the P function varies with r is illustrated in figure 10 for 
a simple one-dimensional design space. The single design variable is x, the objective 
function is f(x), and the single inequality constraint is gj = x - b > 0. The minimum 
value of P for r = r j is Mj; the minimum value of P for r = ^ is M 2 , and so 
forth. Figure 10 is based upon figure 4(a) in reference 12. 


For each value of r an unconstrained minimization is carried out. This procedure 
involves several two-step searches. First, it is necessary to generate a set of direction 
cosines Z which defines a move direction in design-variable space. Second, a one- 
dimensional search is carried out to find the minimum value of P in the Z direction. 
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Figure 10.- Variation of P function with x for several 
values of the constraint weighting factor r. The mini- 
mum value of P/r^ denoted by M^. This is based 

upon figure 4(a) in reference 12. 

In the present study, the direction-generating algorithm which was used almost 
exclusively was the Davidon- Fletcher- Powell method. Once a set of direction cosines Z 
has been found, the design variables C can be expressed as 

C = C* + S • Z (B25) 

in which S is the distance traveled in the Z direction and C* is the vector of design 
variables for S = 0. The search in the Z direction is called a one-dimensional search 
because S is the only unknown. 

In the one-dimensional search used to obtain the present results, the constraints 
gj(c ) are approximated by polynomials in S. The polynomials never exceed second 

degree. Such a procedure is described in references 19 and 20. The flutter speed was 
found to be especially well- suited for polynomial approximation. 

A more complete discussion of these optimization techniques is contained in refer- 
ences 12 and 17 to 23. 
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DESCRIPTION OF COMPUTER PROGRAM 

The computer program SWIFT is written in FORTRAN IV for the CDC 6000 series 
computers at the Langley Research Center. The program also makes use of CDC library 
subroutine MATRIX (ref. 24). For these reasons the unmodified program may not be 
operational on other computer systems. A listing of the computer program is included 
so that mathematical operations can be examined in detail and so that modified programs 
can be developed. In the, latter case this appendix can be considered to be an abbreviated 
user’s document. It is assumed that the user is familiar with the material contained in 
the preceding sections of this report. A flow chart of the program is presented in 
figure 11. 

SWIFT can be used to analyze or design a cantilever plate-like wing. The program 
calculates stress and deflections, natural frequencies, and flutter speeds and frequencies. 
When used in the design mode, the program can be used to find the minimum-weight wing 
cover-thickness distribution which satisfies flutter, strength, and minimum-gage design 
requirements. A six-design variable -strength design requires 4 to 6 minutes central- 
processing-unit execution (CPU) time. A six-design variable -flutter design requires 5 
to 7 minutes CPU time. A strength-flutter (combined) design requires 6 to 8 minutes 
CPU time. 


Description of Input 

All of the input data are described in the main program. Most of the data are in 
namelist form. Additional description of input is given as follows: 

The first card contains an optional word that is printed following the words 'WING 
MATERIAL" in the output. "TITANIUM" is used in the sample problem in this appendix. 

IOPT=l This option is used to find a value of OMEGA, the flutter fre- 

quency, and V, the flutter speed, in the neighborhood of the flut- 
ter solution. Real and imaginary values of the complex determi- 
nant are calculated for values of OMEGA and V starting with 
OMINIT and VINIT and ending with OMFIN and VFIN in increments 
of DELOM and DELV. A solution occurs when both the real and 
imaginary parts of the flutter determinant are zero simultaneously. 
Natural frequencies can be calculated by setting RHOA, the air 
density, equal to zero. The effect of aerodynamic damping on the 
vibration frequencies can be explored by setting only V equal to 
zero. 
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GEl> 


Call FLUTTER & 
STRESS and exhibit 
matrix construction 

Call aUTTER 
and print grid 
of OMEGA and V 


Read: 

WMAT 
$ OPTION 
$ WING 
$ AIR 
$ CASE 
$ FLUTER 



Call PENFUN, FLUTTER, STRESS, 

& MINGAGE to determine constraint 
values for initial design 


Call FPOWEL to calculate 
new move direction 


Call ONEDIM to perform 
one-dimensional search 



Figure 11.- Simplified flow chart of computer program SWIFT. 
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IOPT=2 

IPRINT=0 

IPRINT=1,2, or 3 
ITfflCK=2 

DC (I) 

NC 

IHP=0 

IHP=1 

NN 

NCS 

DELOM, DELV 


Once the flutter-solution neighborhood has been established with 
IOPT=l, the program can be rerun with IOPT=2 to calculate flut- 
ter speed and frequency very accurately (within the limits of the 
theory) and to perform an optimization. 

This option shows a printout of matrix construction. Matrix print 
statements require that NN, the number of spanwise stations, be 
set equal to 4. Final results can be used in deflection analysis. 

These print options are used in the optimization mode. 

The ITHICK=2 option is considered to give a more practical thick- 
ness distribution than the ITHICK=1 option. It is suggested that 
C(l) be set equal to the minimum-gage requirement. Limit active 
design variables to C(10), C(ll), C(12), C(14), C(15), and C(16). 
Active design variables are those variables that are adjusted 
during the design process and are identified by nonzero DC(I). 

Set C(2) to C(9) and C(13) equal to zero. 

These quantities are finite difference increments used to calculate 
derivatives with respect to the design variables C(I). If it is 
desired to hold a design variable fixed throughout the design pro- 
cess, set the DC (I) corresponding to that C(I) equal to zero. 

Subscript of the C(I) having the largest subscript. For example, 
if C(l), C(5), C(6), and C(14) are used to describe the thickness, 
then NC = 14. 

Indicates that wing-depth-variation effects (sometimes called 
thickness effects) are not accounted for in piston theory. 

Indicates that wing-depth-variation effects are accounted for in 
piston theory. 

Number of semispan stations, not counting the root station, where 
10 is maximum and 6 to 8 appear to be adequate. 

Number of chordwise stations used in numerical integration 
divided by 10. For example, when NCS=1, 10 stations are used 
by subroutine MGAUSS. 

Increments for OMEGA and V. Although these increments can 
be large for IOPT=l, they must be small for IOPT=2 where they 
are used as finite difference increments in a Newton-Raphson 
search to converge on OMEGA and V. For example, DELOM=l 
and DELV=1000 are reasonable numbers for the IOPT=l option. 
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However, for the IOPT=2 option these increments must be 
reduced to the order of DELOM=0.0001 and DELV=0.1. 

R Constraint weighting factor used in penalty -function formulation. 

Choose initial R so that the initial value of the penalty function is 
of the same order of magnitude as the objective function, which is 
the wing weight. In the program, both the stress and minimum- 
gage constraints have a built-in additional multiplier of 30. This 
additional multiplier was inserted because both the stress and 
minimum -gage constraints are averaged over the surface. 

RDC Factor by which R is divided at the end of each unconstrained 

search cycle. Numbers on the order of 5 to 15 are reasonable. 

The larger the value of RDC, the fewer the number of uncon- 
strained search cycles. However, large values of RDC may also 
lead to false indications of convergence. Adjust RDC so that at 
least three unconstrained searches are carried out. 

DESCON Three logical variables. The values of each are T or F indicating 

true or false. The three variables refer to the three constraints 
being considered: flutter, stress, and minimum gage. For 
example, TFT indicates that the user wishes to consider the flut- 
ter and minimum -gage constraints but wishes to ignore the stress 
constraint. When a constraint is ignored, no analysis is performed 
to check that constraint. 

Operational Hints 

Flutter .- There are many values of OMEGA and V which satisfy the equations 
obtained by setting the flutter determinant equal to zero. During a computer run in the 
optimization mode IOPT=2 , the program locks in on a flutter solution and follows it 
throughout the optimization process. To obtain the initial solution — that is, OMEN and 
VIN — the user runs the IOPT=l option and chooses the solution which has the lowest flut- 
ter speed. For the final design, however, the program may not be examining the flutter 
solution corresponding to the lowest flutter speed. 

The only way to guarantee that the final design does satisfy the flutter constraint is 
to rerun the IOPT=l option using the final values of C(I). Although experience indicates 
that this procedure is seldom required, it is recommended until a user becomes familiar 
with the particular class of problems on which he is working. 

Convergence .- Convergence to a global extremum is of prime concern to the user. 

In general, the best way to increase one's confidence in a particular solution is to rerun 
the problem with a different initial design, different initial value of R, and/or a different 


32 



APPENDIX C — Continued 


value of RDC. A technique is presented in reference 25 for warping design variable 
space in such a way that the second run will have a greater chance of locating another 
minimum, if one exists. 

The smaller the final values of a constraint function G(I), the closer the design is 
to that constraint. The following values of the constraint functions G(I) are reasonable 
for convergence: G(l) = 0.03, G(2) = 0.15, and G(3) = 0.15. Designs having a value of 
G(I) larger than those indicated should be examined for false convergence. 

Scaling of flutter determinant .- The object of subroutine DEUPPS is to evaluate the 
complex flutter determinant as rapidly as possible. Because of the overriding need for 
speed, it was decided to eliminate certain scaling logic that would normally prevent over- 
flow and underflow. Scaling must be carried out by the user. The elements in the com- 
plex flutter matrix are scaled in statements 101, 102, 103, and 104 in subroutine 
FLUTTER. 

Program listing and sample output . - The remaining part of this appendix contains 
a listing of the program denoted as "SWIFT," a printout of a computer run used to deter- 
mine the flutter -solution neighborhood, and a printout of an optimization run. 


33 



APPENDIX C — Continued 


o 

o o 

o o 

c o 

o 

o 

o c 

o 

o 

o 

o 

o c 

O o 

o o o 

o o o 

o 

o 

o 

o 

o o 

o c 

o o 

o 

o 

o c 

o 

o 

o 

o 

o o 

o o 

c o 

o o 

o o 

o 

o 

LU 

o 

o o 

o o 

o o 

o 

o 

o o 

o 

o 

o 

»H 

fM Cl 

*4* to 

r* 

GO 

O' 

O 3 

<M 

m 








*—4 r-4 

r—i 

H 



OOOOOOOOOOOOOO 

OOOOOOOOOOOOOO 

OOOOOOOOOOOOOO 

ococoocooooooo 

OOCCOOOOOOOOGO 


OOOOOOOOOOOOOO 

OOOOOOOOOOOOOO 

OOOOOOOOOOOOOO 

OOOOOOOOOOOOOO 

OOOOOOOOOOOOOO 

OO^f\Jrn-4-ir\'0h-crCT'O^rg 




3 

LL 


M 
2 
»— * 




































uj o 

>■ 

X 
































X 




t/> 


CD 
























« 









tn 





O 


CL 
































* 




a- 

2 

Q 

n 



V 




















• 









* 




co 

*-i 

UJ 

3 



if 




















2 









< 




UJ 

-J 

*— « 




• 


































-J 

3 

2 


V. 

^0 




•> 
















W 









— 




CD 

< 


O 


* 

in 




















X 









nS- 




< o 

o 

V— 1 


• 

o 




X 
















— » 













*— 1 


UJ 

3 


< 

»—4 




X 
















CM 








X 

w 




QC 

UJ 

Q. 

3 


1-4 

a 




X 
















# 








« 

0 




< 

X 

CO 

CD 


CL 

UJ 




« 
























-s. 


X 



>>- 


»— < 


UJ 

2 




X 




X 


X 










• 








to 

4- 

CD 





CO 

ac 


3 

3 




0 






CM 










O 








3 





3 

UJ 

< 

3 


M 

“5 




00 






** 










3 








X 

(M 

— 



3 

3 


CO 


QC 





» 






•> 






X 




CO 








a 

# 


H- 


a. 

< 

UJ 

1-4 


O 





* 




CM 




X 


X 


f- 




* 



X 





2 

* 

* 

3 


2 

a-4 

o 

o 







2 




* 


(O 


in 








e 



0 





*— 

< 

# 

a 


*— » 

3 

2 



X 

♦ 

s 



0 

* 



* 


UJ 




>r 


a 




CD 



-4- 




* 


* 

< 

2 



►—4 

UJ 

CO 


3 

CO 

V 



►— 4 

2 



• 


UJ 


a 


a 


CD 




3 



a 




X 


— 


3 


UJ 

2 

3 

CO 


O 

QC 

•‘v 



K* 

a 



2 


QC 


* 




• 







* 




X 


m 

— * 

II 


X 

*-« 

O 

UJ 


2 

UJ 




CL 




1-4 


O 


• 


• 


vO 







a-4 




X 



CO 

1 A 


3 


UJ 

2 


UJ 

o 

• 



*—< 

h- 

X 




UJ 


2 


2 


r4 




X 



CO 


X 






3 



Q 

CO 

x 


QC 

2 

O' 



a: 

*— • 

CM 


X 


a 


»— * 




UJ 




m 



Q. 


cn 



oc 

0 


a 


3 

2 


o 


3 

< 

o 



0 

O 

lf^ 


h- 












CM 








< 


0 

< 


O 

< 

UJ 

1—4 


co 

CO 

O' 



CO 

z 

» 




* 


* 


* 


II 



X 

» 



* 


a 



3 

4 - 


►*- 




3 

X 



1 

H 



UJ 

0 CD 


00 


CD 


OD 


00 




X 

* 

CD 



00 


00 


• y~ 


4- 

» 


3 

a 

< 

3 


O 

2 

» 



0 

O 

• 


• 


0 


• 


• 


oc 



9 

9 



• 


• 


0 


<1 


H 


CO 

CO 

»— 4 



2 

*-4 

> 




* 

'O 


>0 


>0 




>c 


Q 


in 

Z 4) 



sO 


O 


3 

1 

♦ 

CM 

o 


O 3 

OC 

LL 


< 

UJ 

<r 



O 

X 

H 


H 


*~4 


»— ♦ 


f~4 


• 


a 

*—4 

4—4 



3 


3 


a 


•M. 

if 

(N 


X 

2 

a 

a 



H 

2 : 



2 

O 

UJ 


UJ 


UJ 


UJ 


UJ 


• 


O 


3 



3 


3 



to 

CM 

if 

II 



< 

a 



OC 

c/) 




»-4 

« 


* 




♦ 


♦ 


• 


3 

* 

* 



* 


* 



<0 


< 

►— 


3 

3 

QC 

2 


UJ 





2 

•» 

M 


H 


M 


II 


II 


• 


< 

H 

II 



II 


II 


0 

3 

w 

if 

3 


2 

CO 

a. 

O 


3 






* 











• 


* 

• 








9 

2 

0 


CL 


a-4 

2 

a 

4—4 


3 





* 

X 

CL 


<r 


O 


CO 


0 


9 


• 

40 

2 



3 


3 


co 



r^ 

3 


oC 

a 

< 

3 


3 

<► 




>s. 

X 

< 


UJ 


2 


<0 


QC 


• 


O 

3 O 



• 


X 


3 

O 

♦ 


3 


CL 

u 


< 


-J 

O 

CO 



* 

X 

• 


QC 


< 


• 


• 


• 


• 

* 

X 



• 





3 



O 




2 

2 


LL 

>r 

CO 



X 

e 



< 


• 


• 


• 


• 


a 

If 

CL 



• 




• 

X 


0 

•> 


o u_ 

<1 

4—4 



nO 




X 

* 

• 




• 


• 


• 


• 


• 


9 



9 




* 

3 

H 


3 


2 

Q 


X 


Q 

m 




X 

X 

• 




• 


• 


• 


• 


0 

O 

• 



• 




X 



4- 

o 

• 

< 


2 

QC 


UJ 

1 

O 



# 

CM 

• 




• 


• 


• 


• 


• 

QC 

• 



• 




0 

QC 

3 


CM 

x 


QC 

♦H 

UJ 


ao 

2 

QC 



X 

O 

• 




• 


• 


• 


• 


• 

O 

9 



• 




O' 

O 

O 

C 

H 

< 

o 

UJ 

• 

3 


>— * 

V- 




m 

H 

O 




• 


a 


• 


O 


• 

X 

• 



• 




X 

3 

X 

if 

3 

oc 

<r 

CO 

to CO 

UJ 


QC 





0 

X 

© 




a 


• 


• 


»— 4 


0 

O 

• 



• 




* 



— 


O IL! 

X 

uj Z 

O 


O 





r-H 

* 

• 




UJ 


• 


• 


h» 


• 


• 



• 




X 

2 

•4- sO 

a 

o 

oc 

3 

2 O 

if 


to 

o 





> 

9 




3 


• 


• 


< 


0 

y- 

O 



0 




X 

O 

X 


z 

QC 


2 

*— • M-4 

X 


UJ 

< 





QC 

0 




O 


• 


• 


QC 


• 

0 

• 



1/) 


C 

X 

X 

3 

ro 


*-4 

a 

c 


3 3 

in 

if 

QC 

2 

Q 



X 

< 

• 




2 


• 


• 




3 

0 

• 



3 

X 

3 

-* 


H 

■If 

O 

3 


3 

< 

3 a. 

X 

O CL 


3 



X 

O 

O 


< 


< 


• 


• 


X 


< 

QC 

> 


X 

3 

* 

h 

X 


<T 

II 


> 

2 



O Cj 


2 


1 

a 



X 

2 

t— • 


UJ 




• 


0 


1— 


•—4 


h~- 


# 

3 

X 

< 

X 


3 


4- 

o 

M 

CO 

OC 

QC 


M 

O 


a 



# 

3 

H* 


QC 


a 


• 


Q 


CL 

X 

OC 


a^4 


X 

Q 

X 

a. 

X 


O 

QC 


CD 

< 


o 

CD K 


2 

2 

LU 

»- 

— 


X 

O 

<1 


< 


LU 




QC 

«— 

3 

* 

3 

X 

CO 


X 

C 

X 




3 

<- 

— > 

o 

x 


U. 

3 3 



*— 4 

o 

CO 

O 

— > 

m 

CD 

OC 




UJ 


Z 

X 

O 

*■ 

O 

x y~ 

X 

2 


X 

X 


CD 




CD 


3 


2 


CO CL 


3 

> 

z 

« 

H 

0 

c 




UJ 


2 

X 

<L 

* 

X 

X 

* 

X 

< 

X 

Ul 





2 



# 

♦- 


CD 

LU 

o 

CO 

2 

X 

X 

LL 

UJ 

QC 

<1 

f4 

w-i 

O 

l— 


0 X 

CO 

* 

CL 

X 0 

X 

X 

X 

X 


0 



to 


a 



X 


w* 


X 

»— i 

LL* 

QC i-4 


o 

CO 

CL 

UJ 

•fr 

< 


LU 

LJ 

X 

<r 

*■ 


X 

CO 

X 


X 

in 







C5 


CO 



X 

■If 

0 

X 

3 

3 

3 

UJ 



M 

UJ 

h* 

UJ 

w 


X 

UJ 

■fr 

LL 

X 

• 

X 

•—4 

X 





O 


to 



2 


CO 



3 

3 

3 

<1 


o 

_) 

X UJ 

h* 

UJ 

y~ 

LL 

X 

V— 

v— 


k— 4 

Q. 

X 

OC 

X 

Ul 

X 

2 : 


c 


y- 


2 


CO 



3 





V- 

V— 

+ 

a 

CO 

2 

< 

3 x 

<X 

2 

< 

UJ 

LU 

< 

< 

<1 

LL 

CO 

X 

3 

X 

• 


UJ 


c 


< 


Mu 4 


< 



O 


u 


• 

< 

<1 

X 

O 

*— i 

X 

> 

C 3 

X 


CO 

OC 

O 

0 

X 

X 

* 

< 

X 

CO 


3 


CO 


CL 


2 : 


2 


X 



> 


Q- 


* 

X 

X 

Nt) 

o 


LL 



QC 

X 

if 

* 

* 

* 

CL 

oc 

X 

* 


* 


* 


* 


# 


QC 


•* 


* 

X 


» 


* 


X 

QC 

oc 


oc 

l/) 




o 

3 

X 

X 

X 

X 

O 

0 

O 

X 

0 

X 

© 

X 

c 

X 

© 

X 

© 

O 

0 

X 

0 

X 

<M 

c 

X 

c 

X 

0 

3 

O 

0 

X 

a. 


CO 



LL 

X 

in 

in 

in 

in 

LL 

LL 

f— 4 

in 

* 

in 

* 

in 


in 

* 

m 

* 

LL 


m 

■fr 

in 

CM 

* 

u> 

* 

iri 

* 

O' 

3 

3 

ao 


I 

3 




H 


m 

<£r 

in 




r4 

r~4 

CM 

CM 

co 

CO 



in 

in 


vO 1^* 


oc 

CO 

co 

O' 

O' 

3 

3 

CM 



H 


3 

*— • 



3 






O 

cm 












CM 












CO 

3 













O 

O 












3 












O 

<n 













^4 

IT 












in 












in. 

0 

in 



o 

o o 

c 

c 



































34 



APPENDIX C — Continued 


o o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

O 

coo 

o o 

o 

o 

o 

o o 

o o 

o o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o © 

o o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

O 

o o o 

o o o 

o 

© 

o o 

o o 

o o 

o 

o 

o 

o 

o 

o 

o 

o o 

© 

o o 

o o 

o 

o 

o 

o 

o 

o 

o 

o c 

o 

O 

o o o 

o O 

o 

o 

o 

o c 

c o 

o c 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o o 

O O 

o 

o 

o 

o 

o 

o 

o 

o o 

o o 

o o o 

o o 

o 

o 

o 

© o 

o o 

o o 

© 

o 

o 

o 

o 

o 

o 

o o 

Q 

o o 

o o 

o 

o o 

o 

o 

o 

o o o o 

o 

o o o 

o o 

o o 

o 

o o o o 

o o 

o 

o 

© 

o 

o 

o 

© 

o o 

O 

o o 

m <f 

tr> 

o r- 

00 

CP 

o 

H 

eg re 

m 


CP O *«4 

eg 

ce 

in wo h- ao <p 

o 

^4 

eg 

re 


in 


r*- 00 

CP 

o 

< 


4- 



m 

in 

m in 

IT. IT, 

tri m 

in o 

'O o 

o 

o 

>o o 

sO o 

r- 

r- 

p- 


h- 

r- 

r- 

r- h- 

h- 

00 00 














• 






* 

















* 















NO 


rr» 




• 


r> 















X 


eg 










*> 



r4 






o 


— 















o 






V. 






* 



ID 


eg 




LU 


>*- 















oc 


r> 




m 






CO 



* 






CO 


* 















<1 


H— 




* 



Vs 



h- 



M 


W 




V. 


* 















UJ 


3 




* 



Vs 








s^. 




• 


• 










V 





oo 


n 







CD 



UJ 





O 




I'- 


z 








Vs 

Vs 

V. 







i— 




* 


Vs 

* 



X 



# 


< 




LL 


*— • 








* 

* 

•tt- 





X 


z 




«•* 


— • 




UJ 



eg 


* 






w 








0 

0 

9 





UJ 


r-4 

V. 



st 


CD 

CD 



oc 



*—* 


X 




* 


Vs 





V 

«*. 


o 

o 

o 





-J 


«: 

V V 






* 



o 



* 


o 




00 







re 

re 


UJ 

UJ 

LL 

•v 

V 

— 


CL 


a 

eg re 



«<-r 


CD 

CO 



z 



w 


r-4 




• 


eg 





— i 

t—4 


CO 

CO 

CO 


# 



X 



# * 



o 


1 

1 



*— < 



Q 


• 




o 


# 





# 

H 


Vs 

V, 

Vs 

0 

0 

• 


u 


h- 

* # 




Vs 

r-H 

r-4 






< 


OO 




r-4 







ii 

II 


00 

00 

OO 

(_) 

a 

o 


o 


oc 

CO cc 

v* 


+ 

CO 


-«r 

Vs 


UJ 



* 


• 




LU 


# 








z 

z 

z 

UJ 

Ul 

LU 




o 

* * 




* 



Vs 


o 



X 


40 




* 


o 





z 

CO 


< 

< 

< 

00 

CO 

00 


1 


X 

— . — • 

— 


eg 

— 


— * 

— * 




in 

o 


H 




II 


UJ 





z 

a 


M 

»-* 

*— • 

Vs 

'v 

V. 




CO 

re rO 

UJ 


* 

m 

m 

ro 

UJ 


UJ 


UJ 

r-4 


ID 






CO 






z 


o 

c 

c 

0 

• 

0 


eg 



* * 

h- 



« 

* 


h— 


oi 





« 




• 

Vs 

* 






# 


< 

< 

< 

Z 

Z 

Z 




1 

* * 

♦ 


< 

tf 



4* 


UJ 


CD 

00 


II 




O 

> 

• 


X 




# 

Vs 

oc 

oc 

CC 

r-4 

r-4 

t— * 


# 



< < 

Ol 



<t 

<x 

<1 

UJ 


u. 


< 

• 






UJ 


CD 


o 




• 

Vs 








> 


r-4 

* * 

-J 



* 

* 

* 

-J 


u. 


*— < 

o 






00 

CD 

-J 


o 




«•» 

* 

* 

* 

* 

« 

« 

* 





— * — * 

w 


m 


— * 

— 

w 




oc 

r-4 


* 




Vs 

♦ 

*— 


* 




oo 

Z 

ao 

oo 

00 

00 

00 

CD 


Q 


• 

eg no 

Vs 



00 

eg 

vD 

•v 


Q 


< 

UJ 


eg 




0 

< 



0 




o 

O 

• 

• 

• 

9 

0 

0 


Z 


z 

H H 

X 


w 


H 

r-4 

X 




> 

* 


*-r 




z 


X 


I'- 




z 

1-4 


o 

o 

o 

'O 

vO 


<3 


D 

■»»»■ 

* 


o 

nr 

%■» 

— > 

* 


UJ 



M 


* 




M 

II 

re 


LL 





h— 

r-4 

r-4 

•-4 

r-4 

r-4 

r—4 




r-. 

o U 

eg 



o 

o 

u 

eg 


h- 


z 



w 






eg 






X 

C 

LU 

UJ 

UJ 

UJ 

UJ 

UJ 


<1 


t— 


1 


+ 




1 


*— • 


o 



u 




* 

UJ 

» 







a: 

* 

* 

» 

* 

# 

* 


o 


o 

+ + 




+ 

•f 

+ 

rH 


z 


t— « 

* 


Q 




00 

o 

CD 

Vs 

CD 




o 

UJ 

II 

II 

II 

It 

II 

II 


UJ 


3 




< 






n— « 


CO 

eg 


* 




0 

o 

0 

ao 

0 




r-4 

K* 








X 


ac 

eg <nj 

It 


* 

eg 

eg 

eg 

II 


LL 


UJ 

M4 


X 




<£> 

UJ 

>o 

• 

vO 




w 

*— * 

K 

X 

z 

t- 

> 

z 


a 


b- 

* * 




* 


* 



* 


Q 

« 


o 




r-4 


r-4 

NO 

r-4 






►-r 

o 

r-» 

*— « 

-1 

1—4 




OO 

* * 

00 


eg 

* 

* 

* 

CD 


X 


* 



r-4 




UJ 

© 

LU 

^4 

UJ 




II 

U. 

Z 

-J 

u_ 

z 

LU 

LL 


z 


Z 

< < 




< 

< 

< 



h- 


X 

© 


•» 




* 

z 

* 

UJ 

* 





O 

r-4 

UJ 

X 

r-4 

CJ 

> 


u 


a 

# * 



w 

* 

* 

* 



r-4 


fv 

LJ 


OJ 




II 

r-r 

II 

* 

II 




CO 


X 

o 

LJ 

> 






o 




O 



a 



•> 


o 

# 


• 





o 


II 




CO 

Z 

b- 

O 







u? 



m 

o 



r- 

r-4 

in 

o 


# 


i-4 

X 





Vs 

< 

<T 

< 


• 



z 

O 

Z 








t- 


X 

H ^4 

z 


+ 


H 

r-4 

Z 


CO 


*V 

o 


r-4 



* 


LU 

o 

< 

• 



a 

*-4 

LU 








< 


r-4 

•— w> 

< 



nr 

«— 

* 

< 


UJ 


# 

r-4 


UJ 


X 

X 


-J 

X 


• 



►H 


X 







Vs 

CL 


ac 

© o 



—> ■ 

u 

u 

o 



_J 


UJ 

r, 


# 



X 



CL 


0 



1- 

<r 

UJ 








LU 


i— 




rJ 






CD 

00 

> 

on 


II 


r- 

X 


X 



UJ 



< 

t~ 

oc 







* 

b~ 


<1 

+ + 




4- 

■4* 




< 

UJ 


• 




it 

* 


r-4 



a 



t— 

00 

o 







CO 

1-4 


X 


00 


w 




00 


r-4 

-> H 

'O 


— 


oo 

X 


eg 



z> 



co 


z 







z 




< < 

oo 


o 

< 

< 

< 

00 


QC 

M 

O 

r-4 




uj <n 

0 

•> 



h- 


V. 


UJ 








o 

1 


1 

* * 

Vs 


X 

« 

* 

♦ 

Vs 


< 

K 

<r 

UJ 


eg 



o 







V Z 

oo 








r-4 




—s A 

>- 


rH 




>* 


> 

< 


* 


M 


K 

,-H 


X 


X 

K 


* 

< 

1— < 

Q 




0 



t— 

r— 1 


o 

o 

1 

— * 

* 

nO 

o 

4- 

1 



> 

CO 

II 


* 


or 

Vs 


X 


X 

-J 


in cl 


z 







a 

* 


* 

H rH 

rH 

4P, 

X 


H 

r-4 

r-4 


Z* 

M 

uu 



w 


UJ 

# 


X 


X 

< 


z oo 

a 

< 


h- 





a 

X 


X 



H 

re 



w 


eg o 

QC 

b- 



o 


a 

X 

O 






O 

9-4 

oc 



z 






in 


m 

o o 

II 


* 

u 

o 

o 

II 


r-4 

UJ 

< 



« 


o 

X 

z 




O 


r-4 

X 

c 

UJ 

-J 

UJ 





h~ 

» 


» 


II 

It 





M 

O0 

Q O 

eg 


X 


a: 

X 

X 


> 


z 


t- UJ 

X 

o 

< 

X 





3 

eg 


eg 

+ + 

< 



+ 

+ 

+ 

< 


UJ 


t-H 

M 


eg 


a. 

• o 

II 

K 




< 

00 

o 

z 

r— i 

LU 

-J 


H 

o 

CL 

*—» 


r—4 



X 

0C 





^L 

o 

_J O 



» 




l/) 


»—* 

• o 


K 



< 

t- 

a: 

< 


Z 


z 

* 


* 

— 


o 

< 

A 




O 

* 

< 

z 

— - 


# 


or 

0 


z 

in 

z 

z 


00 

LL 

LL 

oc 

M 

o 

z 

-J 

Cu 



II 


tl 

O' fp 

UJ 

r-4 

OJ 

in 

CP 

cn 

UJ 

M 

X 

© 


u 


X 


*—* 


LL 


z 

9—4 

o 

— 


O 

o 


z 

z 

•—< 

c 

X 






r«4 

C* 

X 

K 



r-4 

oc 

X 

in 



It 


r4 


< 

• o 

Q. 

UJ 


Cl 

ao 

* 



* 

M 

r-* 

LL 

r-4 

LU 

-J 

-If 

b- 


t— 

w w 

LU 

K 


nr 

w 

<« 

UJ 

h- 

V. 

K 


X 


X 





CO 

Q 

«• 

00 

9 

Vs 

OC 

a: 

Vs 




y- 

oc 

<1 

Vs 

CL 


z 

u o 

X 

r-4 

* 

o 

o 

u 

X 

M 

Vs 

> 

» 

in 


eg 


■fr 

• 

o 

1— « 


» 

LL? 

o 

Vs 

ID 

LU 

V 

< 

< 

< 

9-4 

o 

Z 

V 

O 


r-« 

— 

3 

w 

** 

w 

w 



W 

— r 

-J 

# 



W 

CO 

*— 


UJ 

X 

to 

r4 

oc 

r4 

W. 

CD 

CD 

W 

o 

© 

© 

z 

z 

►— 

w 

r— t 


OC 

♦ + 

X 

X 


+ 

+ 

+ 

X 

X 

b- 

< 

X 

K 


H* 

s 

t— 

e 

UJ 

UJ 

00 

♦ CL 

UJ 

b~ 

X 

X 

h- 

LU 

LL 

UJ 

S-4 

r-4 

LL 

1— 



CL 

X X 

re 

eg 

< 

X 

X 

X 


eg 

<1 

Z 

r-4 

< 


< 

4D 

< 


Q. 

00 

< 

IT. 

LJ 

* 

< 

3 

3 

<1 

X 

X 

X 




<a 



r-4 

O' O' 


H 

X 

*o 

ej 



r-4 

X 

< 

X 

X 


X 

r-4 

X 

# 

CO 

* 

X 

LL 

U 

II 

X 

Z 

Z 

X 

o 

c 

o 

> 

> 

> 

X 




ss- *■ 

X 

X 

cr 


IT 

m 

X 

X 

CL 

* 

r-4 

or 


0£ 

UJ 

CL 


# 

X 

# 

X 

«■ 


CC 

* 

* 

oc 


* 


* 

■fr 

■» 

ac 



* 

X X 

eg 

vO 

O 

X 

X 

X 

eg >c 

C 

X 

r-4 

c 

— * 

C 

# 

O 

# 

X 

in 

X 

CO 

X 

CD 

c 

X 

X 

a 

X 

X 

X 

X 

X 

X 

O 

X 


X 

00 GO 

rW 

eg 

Li. 

00 

00 

cc 

r-4 

eg 

LL 



LL 

CO 

LL 

II 

LL 

■* 

in 


m 

r-4 

in 

* 

LL 

in 

in 

u. 

u- 

in 

in 

in 

IT. 

tn 

LL 

in 

Vs 

IT. 

r-4 H 

r-4 


eg 

r-4 

r-4 

rH 

r-4 

r-4 

H 

r-4 

eg 

CO 

r-4 


r4 

>i* 

r-4 

r-4 

eg 

eg 

eg 

re 

re 

in 

r-4 

eg 

eg 

r-4 

eg 

re 

<*■ 

in 

vo 

O 

•-4 

^4 

r4 




m 






m 



m 


re 


o 











eg 







in 







O 






m 



u> 


in 


u\ 























35 


1- INTERMEDIATE PRINTOUT, 3 - EXTENDED PRINTOUT*/ 8200000 
75 X* I SN OWCH =*I 2, 5X*0 - SCL ID WING, 1 - SANDWICH CONSTRUCTION*/ 8300000 
85X* IHP =*I 2* 5X *0 - THICKNESS VARIATION EFFECTS NOT ACCOUNTED F 8400000 
8 OR IN PISTON THEORY, 1 - THICKNESS EFFECTS ARE USED*/ 8500000 
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SUM( NN )=R< JJ)*FOFX(NN)+SUM<NN) 31100000 
00 7 I JK= 1, NUMBER 312CC000 
SUM ( I J K ) = SUM( I JK )/DELTA 31300000 
RETURN 31A0C000 
FN0 31500000 



SUBROUTINE ANDER 31600000 

31700000 

ANALYTICAL DERIVATIVES OF WEIGHT WITH RESPECT TO DESIGN VARIABLES 31800000 


APPENDIX C — Continued 
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ADC 121 = A0C4I - ADC 12 ) 34600000 
RETURN 34700000 
END 3480000C 
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IF(KCERIV.EQ.O) GO TO 50CC 18270C000 

0MEGA1 = OHSAVE 1 828CC000 

VI = V SAVE 182900000 

5000 RETURN 183C00000 

END 183100000 
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SUBROUTINE SMAT 189200000 

189300000 

GENERATES DIAGONAL ELEMENTS IN S MATRIX 18940C000 

189500000 

REAL MU 189600000 
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END 199200000 



APPENDIX C — Continued 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CJ 

o 

o 

a 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

o 

CJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

u 

o 

o 

o 

CJ 

Q 

o 

o 

V 

u 

o 

o 

o 

o 

o 

o 

o 

o 

UJ 

o 

o 

o 

o 

o 

CJ 

o 

o 

© 

© 

o 

a 

o 

o 

a 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CO 

4* 

in 

4* 

r* 

03 

O' 

o 

H 

CM 

m 

4- 

in 

o 

r- 

cc 

O' 

o 

rH 

CM 

ro 

4- 

in 

vO 

h- 

00 

O' 

o 

rH 

CM 

co 

sT 

in 

o 

r* 

00 

O' 

o 

H 

CM 

CO 

N* 

O' 



O' 

O' 

O' 

O' 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

«M 

rH 

H 

rH 

rH 

rH 

H 

rH 

rH 

rH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

co 

CO 

m 

CO 

CO 


ur 

cr 

O' 

O' 

O' 

O' 

o 

o 

o 

o 

o 

o 

o 

u 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


H 

H 

H 

rH 

rH 

rH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

(M 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

«M 

CM 

CM 

CM 

CM 

CM 

CM 

CM 


h- cx 



2 <1 



UJ 



i cl 



O X 



X 

\ 


<r o 

\ 


2 

M 


X < 

00 



• 


(X u. 

M3 



rH 


< 

UJ 



* 

— * 

>- 

II 

H 

<X (X 


H 

< n 

:> 

» 

UJ 


in 

X x 

a 

w 

UJ H- 

2 

a 

X 

< 

X • 

a 2 


r 

X a 


— X 

o v- 

* 

»H ^ 

O CO 

» 

rH x 

1— I 

ao 

♦ 2 

a 

• 

in * 

* 

M> 

CO 

► a 

fH 

LL O 

CO 2 

UJ 

♦ 2 

•-» < 

* 

. r 

* 

II 

rH CO 

X ^ 


rH o 

</> 

< 

♦ a. 



> 

rH 



-H 


— * 




< 

M 



•> 

u 




cm 

X 

ro 



nr 

tx 

HT 



CL 


a 



X 

• 

X 



» 


* 



<r 

4 

< 


MH 

— * o 

— »• 

a 

— » 

* 

— . UJ 

M 

UJ 

— 

CM 


• 

X 



•• o 

CM 

o 

r-4 

It 

—1 * 

Hr 

* 

• 



il 

4* 

CM 

> 

UL rH 

* 

rH 

Hr 


* X 

> 

3 

CL 


> < 


< 

X 


* o 

<[ 

o 

* 

< 

< » 

* 

* 

< 

o 

* > 

<x 

> 


X 

< * 

o 

* 

4* 

CX 

o < 

X 

< 

rH 

* 

X o 

cc 

o 

X 

• 

<X X 

* 

I 


CM 

* CX 

• 

% 

o 

4* 

• * 

4- 


— * 

CM • 

4 

e 

> 

*— 1 

4- cm 

«*H 

CM 

— 4- 


— 4 

M 

4 


CO 

rr «*h 

♦ 

— * 

r- * •—» 

■^r 



* 

UJ 




o 




in 

•> 


to 












— 

rn 

r 

rn 


CO 


rn 


•» 

•» 


U) 


► 

u 




(JU 





n) 


CL 












— 

•* 

CM 

r 


Hr 


r 


CM 

CM 


* 


CO 

1— « 




X 




CO 

CO 


UJ 












rH 

r< 

Hr 

CM 


CO 


<o 

— ■ 

«r 

HT 


CM 


M 

cx 




o 




•» 

2 


•» 












r 

w 

CO 

nr 


<4 


w 

O CO 

UL 


X 


X 

h- 









2 


CO 












rH UL 

4 

UL 

C? <M 


UL 

< 

* 

w 

o a 


in 

< 




l/) 




ni 

r 


CO 












w 


CM 

* 

<1 

X 


* 

X 

CM 

* 

< 

1 


•» 

X 




UJ 




O) 

CO 


•» 












to 

<1 

— X 

<1 

X 

u 


<t 

< 

X 

c 

X 








CO 




» 

CM 




X 










4 

O C3U 

e? 

<t 

1 


o 

r 

o 

o 

< 

* 


X 

to 




X 




CO 

2 




u 










NUJ< 1 

UJ 

r 

Hr 


UJ 

X 

1 

UJ 

•* 

H 


►H 

— 









2 


rH 












X 

X 

X w 

X 

X 

* 


X 

< 


X 

X 

00 


CX 





> 





» 


w 


Q 










a 

O <f * 

a 

<t 

>* 


a 

UJ 

* 

a 

< 

a 


H- 

to 




cc 




< 

CM 


o 


2 










1 

* 

•* <h 

* 

UJ 

oo 


* 

cx ♦ 

UJ UJ 


< 

X 




o 




< 

CM 


« 


00 











< 

X 00 

< 

cx 

a 


< 

< 

00 

< 

oc 

4 


X 

cx 




UJ 





2 


X 


•H 










* 

* 

< a 

* 

< 

UU 


* 

*HT 

a 

* 

< 

-H 



UJ 




X 




A 

2 


X 


r 










4- 

<r 

UU UJ 

< 

w 

4 


<1 

X 

UJ 

< 

HT 

rH 

o 

o 

H 




h- 




co 

#• 


r 


K- 




• 






00 

C CX 4 

o 

X 

rH 

•— 

c 

X 

4 

o 

X 

2 

c 

< 



>r 






CO 

CM 


UJ 

G. 

2 




>*■ 


— » 




a 

X 

< 

X 

X 

CM 

rr 

X 

CL 

rn 

X 

X 

2 


< 

< 


• 

■v 

z 




» 

2 


•> 

X 

►H 




^ ro 

CO 





UJ 

cx 

w »H 

cx 

a 

CM 

n 

cx 

X 

tm* 

cx 

a 

*H4 

►*- 


*-*< 


H 

N 

c 




ro 

2 




cx 




— CO 

CO 

» 

rH 



4 

4 

>r 2 

tf 

X 

2 

r 

* 

o 

r 

* 

X 

r 



h~ 


rH 


1- 




CO 

» 

> a 

♦ 

a 




• 2 

2 

M 

2 



r 

# 

X 2 

• 

o 

2 

CM 

• 

II 

rH 

• 

u 

rH 

a 

V 

(X 


UJ 


00 





rH 


X 

< 

>— i 

CM 



rH 2 

2 

W' 

2 


CM 

M 

CM 

Q_ rn 

CM 

II 

rn 

w 

CM 

r-- 

2 

CM 

II 

2 



UJ 


o 

, 'v 





u 

2 

<t 

cx 

a 

•> 

* 

'T 


4- * 

•> 

< 

r 


CM 

» 

W 

X » 

Hr 

— - 

» 

a 

Hr 

CM 

2 

Hr 

— » 

2 

UJ 


2 


rH 

rH 

CL 


c 


< 

2 

CJ 

» 

X 

K 

* 



X rH 

rH 

< 

rH 

rH 

2 

M 

* 

u — • 

* 

rH 

*— * 

X 

* 

CM 

M 

* 

rr 

M 

2 

X 

#-H 

<— 

X 

X 





<1 

r 

UJ o 

cx 

a. 




< II 

II 

II 

II 

2 

2 

w 


II — 

<4 

2 

rr 

* 

>* 

2 

Hf 

•4- 

r 

Hr 


rH 


>v 

rH 

H 

4 




V 

2 

X 

2 

• 

o 

O CO 


o *H 

•9 

• 

►H 

2 

4 

<r 

CO 

— <1 

oo 

2 

<r 

sT oo 

2 


oo 

rH 

<r 

h 


I ^ 

^r 

^r» 

W 

X 

X 


X 

*- 

a 

<1 

<1 

M 

UJ 

cx 

'O* 

W 




4 

•H 

< 

CL 

rn 

CL 

rn 

<1 

rH 

a 

rn 

<x 

a. 

2 

<1 

X 

h- 

K- 

* 

K 

K 

h- 

X UJ 

2 

X 

V. 


V* 


V 

X 

UJ 

• 

II o 

o 

« 

o 

rH 

II 

It 

UJ 

- II 

UJ 

•> 

It 

X 

UJ 

» 

II 

UJ 

2 

M 

o 

< 

•-H 


<r 

< 

<1 


X 

o 

on 

rH 

CM 

CO 

•4- O U 

M 

rH 

< o 

o 


o 

II 

CM 

X 

II 

rn _J 

II 

rn 

X 

< 

II 

rn 

X 

II 

rn 

X 

CX 

X 

2 

00 

X 

X 

X 


a 

X 

7 

X 

X 

X 

X 

X 

II 

*• 

II 

rH rH 

rH 


CM 

rH 

CM 

< 

<3 

w HJ 

o 


< 

c? 

o 

HT 

< 

o 

-r 

< 

ao 

CX 


> 

CX 

cx 

cx 

<1 

X 

X 

U- 

X 

X 

X 

X 

X 

CM 

00 

X 

X 


u 


2 

2 

UJ 

< 

O UU 

< 

o 

UJ 

* 

< 

o 

ui 

< 

o 

UJ 

X 

o 

Q 

C C 

o 

a 

uj C 

o 

2 

cc 

cc 

cc 

cc 

CD 

X 

CL 

<x < o 

o 

< 

o z 

2 

cc 

X 

< cx 

X 

<1 

CX 

> 

X 

< 

CX 

X 

<1 

cx 

to 

U- 

UJ 

cx 

u_ 

u. 

UL 

CC 

o 

a 

>v 

V* 



"V 

N* 

□ lUOOQ 

o < 

Q 

►H 

k-H 

< 

< 

< <1 

c 

< 


* 

< 

< 

<r 

< 

< 

< 



rH 

CM 







rH 

rH 

CM 

rH 

rH 

rH 
















rH 








rH CM S*- rH 

CM 


o 

o 


88 



l+2.*RH0A*V*V*GAM14*HP(2,n ) 203300000 

AM AG = FPS4* (2«#RH0A*A*0M£GA*F (3 1 1 ) + 2. *R HOA* V*GAM14*0M EGA*HP ( 3 » I J ) 20360C000 

AAC ( INN1,INNU=CMPLX(AREAL,AMAG) 203700000 

AR EAL= AA( INN1, INN22 » +EPS4* (-CM 2* S< 4 , 1 1 +4« #RHOA*A*V*F ( 3« I ) 2C3800000 


APPENDIX C — Continued 


o 

o 

o 

O 

O 

o 

o 

O 

o 

o 

O 

o 

o 

o 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

Q 

o 

o 

O 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

O 

o 

O 

O 

CJ 

O 

o 

O 

O 

a 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

O 

o 

O 

o 

O 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O' 

o 

H 

<M 

m 

4 

in 

<4 

r- 

CO 

o 

o 

f-4 

CM 

CO 

4 * 

in 

o 

r- 

00 

O' 

o 

H 

CM 

m 

m 

4 

4 

4 

4 

4 

4 

4 

4 

4 - 

4 

in 

V\ 

ITS 

in 

in 

in 

m 

in 

m 

in 

o 

o 


sO 

• o 

o 

o 

o 

o 

o 

o 

o 

© 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

(SI 

CM 

CM 

CM 

CM 

(Si 

CM 

CM 

(Si 

CM 

CM 

(si 

CM 

CM 

CM 

CM 

CM 

fM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 














A 



















«— 




A 












»-4 











»— « 












•» 







•> 




•* 












4 







4 




in 












w 







w 




w 












Q. 







a 




CL 












X 





— * 


X 


►M 


X 












* 







* 




# 












< 





• 


<x 


4 


< 












o 





ro 


o 




o 












uu 





w 


LU 


LL 


LU 












X 





UU 


sr 


*»• 


X 












o 



*— « 




o 


> 


o 












♦ 



•> 


> 


* 


* 


* 












4- 



cn 


* 


4 


< 


4 












H 



w 


<1 




♦ 


H 












J 



Q- 


* 


X 


< 


X 












<r 



X 


< 


< 


a 


< 












o 



* 


a 


o 


X 


o 












* 



< 


X 




CL 


* 












> 





cc 


> 


* 


> 












* 



4 


* 




• 














< 



H 


• 


< 


4 


< 








•M 


cn 


a 



X 


CM 


o 


♦ 


a 








fO 


m 


i 



< 


■f 


X 




X 








m 


z 


a: 



o 




cc 




CC 








z 


z 


* 



* 


IH 


•» 


* 


*■ 






m 


z 


» 


• 


<■»* 

> 




• 


in 


• 






m 


» 


H 


<4 


— * 

•f 


4- 


CM 




CM 






z 


H 


II 


+ 





w 


+ 


on 


* 






z 


II 


M 





*-* 


oo 


— ■* 


* 


*** 

o 





► 


►-* 


* 


•-» O 

«k 


* 


*-4 O CM 



< 





H 




— * 


•» 

<1 


m 


CM 


•* <1 

X 


•» 






II 


a 


ro 


4" X 

</) 



X 


4 X 

a 



< 







CM 


CO 




* 

u. o a 


w 

<1 

i 



•* 





•> 


CM 


z 


LL 

•» 

CM 


< 

1 


LL 

► 

w 


LL 

-J 





— * 


Z 


z 


* 

-J 

X 

* 

X 




-1 # 

—» 

«» 

<1 





H 


z 


•> 

— * 

< 

< 

u 

< 

< 

* 

— » 

< 

<x 

4 

— s 

< 

UJ 





z 


•» 


m 

HOW 

1 

o 

•» 

4 

M 

U UJ 

C/5 

M 

CJ? 

CL 





z 


CM 


CM 

•» 

uj a; 


UJ 

-J 

OO 

► 

UJ a 

a 

•* 

UJ 

■4 





•» 


Z 


Z 

f*l £ < 

* 

X < 

a 

m 

X «a 

UJ 

4 

X 

w 





H 


z 


z 

>■» 

u 


4 

O UJ UJ 

■—» 

a 


4* 

W 

a 

X 


z 



II 


II 


H 

a 

* 

X 

00 

# 

a: 


a 

* 

X 

— 

CL 

* 

-J 


CL 



n 


3 



X 

< -I 

Q- 

< 

< 

4— 

X 

< 

CM 

X 

< 

CL 


3 



•» 


* 



* 

* 

a 

UJ 

* 

«•* 

r-H * 

* 

CL 

CM 

» 

* 

X 


►- 





*» 


4— . 

4- < 

X 

+ 

< 

X 

z -t < 

X 

Z 

4 < 

o 


UJ 

m 


"5 




“0 

— i a o 


a _j 

Z C O U Z 

r-c a 

II 


ac m 

> 

* 






X 

ii 

1-4 

x a 

HH 

X 

X 

II 

NH 

X 

X 




z 

» 



M 


H-1 

< 

cc. 


* 

a: 

X 

• 

< 

o c 


• 

< 

CL 

CM 



Z <1 

W 


W 



e> * 

CM 

CM 

* 

o 

rsi o 

* 

H N C 

* 

CM 


o 

(► 

o 

LJ 


o 


o 

* 

• CM 

CM 

• 

II 

CM 

* 

• 

Z (Si 

* 

• 

z 


t m 

LL 

< 


< 


<1 

> 

CM 


z 

CM 

— ► 

z > 

NZZ 

> 

CM 

z 


uj m 

X 

<1 


«x 


<1 

* 

w* 

Z 

z 



z 

* 


** 

z 

* 

— 

*-4 


z 

z 

a 



w 


w 

> 

# 

M 

M 


• 

►— 

> 

* 

* 

»— « 

> 

* 

•* 


0 

z 


w 


W 


W 

* 

4* 

► 

w 

4 (4 

w» 

* 

4 CM 


* 

4 

CM 


>- 


•• 







l/) H 

< 

00 

CM 

< 

< 

to 

CM 

< 

<1 

1/5 

CM 

UJ Z 

•J 




« 


«► 

CJ 

a 

z 

< 

CL 

Z 

<1 

o 

Q. 

Z 


a a. 

Z 

X 


r< 

CM 

CM 4 

CM 

<• 

CM 

X 

UJ z 

It 

UU 


II 

X 

UJ 2 

II 

X UJ 

z 

z 

CL 








CL 

II 


-J 

II 

M 

-J 

OC 

II 

M 

-J 

cc 

II 


►H 

a 

H- 

b- 

h- 

h 

h- 

H- 

K 

» 

LD 

w 

< 

O 

W 

<[ 

* 

o 


<1 

* 

o 

w 

1- 

»— * 

Z 

z 

Z 

Z 

Z 

Z 

Z 

• 

< o 

uu 

< o 

UJ 

t 

< u ui 

• 

< 

o 

z 

w 

t-4 

IM 

t-4 

►— » 

k— i 

•— 


4- X 

< 

cc 

X 

< 

or 

CM 

X 

< 

CL 

4 X 

<1 

C LL 

CL 

CC 

CL 

CL 

CL 

CL 

CL 

♦ 

< < 

< 

<t 

< 

< 

♦ 

< 

< 

< 

•f 

< 

< 

o *-• 

a 

CL 

a. 

a 

a 

CL 

CL 


o 

o 

C\J 


89 



APPENDIX C — Continued 


OOOOOOOOOOO 

OOOOOOOOOOO 

OOOOOOOOOOO 

OOOOOOOOOOO 

OOOOOOOOOOO 

> 3 -in>or~® | 7 'o— 

OOOOOOOOOOO 

(XifMrjfSlrvjtMfMtNifMfMfNJ 


o 


o 

o 


o 


o 

2 T 

« 



O 














UJ 

< 






o 

<N4 






Q 

O 






UJ 

» 





UJ 


CD 

CvJ 

> 



Uj 

o 

iM 

a 

* 



u 

z 

u 


CD 




►-« 

► 

H 

H 



UJ 

o 

< 

o o 



Z 

< 

H 

» 

+ 



»—* 

UJ 

u 

> 


<1 




V 


r-4 

CNJ 


Z> 


Z >0 

m 

u o 

z 

Q 

CO 

a m 




tt. 

a 

Ui 

z *: 

* 

II 

(1 

z> 

cc 

z 

s: _j 

— 1 



»- 

Z> 

M 

Q CD 

CD 

r-» 


u_ 

o') 

LL 

O ^ 

V 

u 

o 

ct 


UJ 

m 

fn 





Q 







O 







90 



APPENDIX C — Continued 


> o o o 
o o o 

Cl O O 

ooo 
O O o 
i-i in cd 

00 00 CD 
OOO 
I CM CM CM 


ooo 
ooo 
ooo 
ooo 
ooo 
cm (d 'J- 
O' O' 
ooo 

CM CM CM 


I o O o o o o 

I O O O o O O I 

O O O o o o 

O O O O O O i 

O O O O O O i 

P- CO O' O H CM l 

O' O' O' O O O i 

O O O r * h( hC I 

CM CM CM CM CM CM I 


OOOOOOOOOOOO 

OOOOOOOOOOOO 

OOOOOOOOOOOO 

OOOOOOOOOOOO 

OOOOOOOOOOOO 

00000>-1*-l*Hi-lr-<r-lr-1 
H H H H H — — f 4 H H H H r~C 

CMCMCMCMCM'MCMCMCMCMCMCM 


X 

a 


H 






ta-t 

QC 


» 






a: 

LL 


in 






H 









< 

a 


a. 


fO 




2l 

Uj 


X 

Cd CM 

m 





K 


•> 

CM CM rr> 

<. 




< 

2 


X — 

CM <t < 






•-* 


<t rH 

CO *H <f 

M 





oc 


21 -H 

CM H H (Cl || II 






a 


2 * 

hC CM < <C II 

«“» 




# 



► in 

< U — — 

fO 




X 

ID 


OO 

r-4 II II — Cd m 

fO 




m 

ac 


O LL 

h( II l| CM cd CM 

z 




* 



Z •* 

<T — — CM Z Z 

2 




m 



— » 

— — fd -c Z Z Z 

• 




*— < 

00 


l/) *H 

ii cm h n z z * • 

m 




* 

UJ 


a «-» 

CM Z Z Z » CM <J- 

m 




X 

ID 


3L * 

— ZZZ • CM CM CO 

2 





»— 4 


•* u\ 

HZ - -(C1IMZZ 

2 




* 

OH 


cn -w 

Z KlIHICZ ZZ 



o 


•* 

*- 


m CO 

ZHTIZZZ-- 

< 


H 


m 

< 


2 * 

*zz:zz-<< 





►h 

2 : 


Z <— 

HZ Z - - < < < 



o 

-o 

X 


H 

•» <n 

w w < «i «i 

O 


h- 

•* 

in 

o 

H 

m ro 

z <* < < <c o o 

K 



F-H 

•» 

z 

•> 

CM I * 

—• -a <t a o i- 



o 

w 

* 


H 

Z O CO 

< O O 1- 


o 

o 

< 

X 

3C 

H 

Z jn 

<10 0 1-1- — — ■ 

fO 



eg 

►h 

o 


► a — 

HI- — Cd <M 

04 



H(\| II 

0 c 

-J 

h 

IMZ< 

q — — CM CM Z 

z 

1- 

in g> 

2 Z 


k- — — • z z z z 

N H CM - Z Z - * 


Z Z — 
+ + < 



r s 

LL 




2 

» 





z 


2 cn 

» 

• 

ro 

m 

CL. 


UJ 




-O -5 



o 



* 


V 

rH 


• 

1- 

re 



Hi Z 

eJ 2 eg 

eg 

eg 

CM 

eg 

CL 


• 

o 




► 


o 


i- 

< o — 

UJ V 



eg 2 




• 

► 

z 

► z 

2 

2 

z 

z 



< 

h- 



ii n 

l-H 


eg 


< 

^ X 

* 

* 


z 

LH 




H 


z 

H 

2 

2 

z 

2 

2 



— 2 





w 


» 



-o — "M 


K 

H 


2 

oc 

ro 






w 





w 

o 


< • 

C Hi 

rH 


< < 

< 


O 

rH 

< 

# CO co v 


<X 

rH 


• CL 

m 



< < «t 

< < 


< 

< 

< 


OC 

o z 

2 


-> < 


o 

2 

< 

'v I • V 

* 

X 



p-4 

M 

w» 



< 

< 

< 

< 

< 

< 

< 

< 

<X 

o 


< o 

z 

2 





eg 

2 


z -o v 

S <1 

c u 


2 

► 

o 












UJ 

rH 

z • 

•> 

» 


— . — 

II 


» 

♦ 

UJ 

H h* H H 

H 


< 


z 

h 

< 

o 











0 

ii 

<r 

*C rH 

rH 


eg cn 



O rH 

z 

x uj x 

i 

o 

z < 


* a 

< 













< 

h- o 

• II 

ii 


• • 

—* 


o 

II 

►— i 

h < in h 

r~» 

2 

o 


2 

o 


II 

H 










2 

2 

< r cj 

o •- 

-o 


o a 

<1 

UJ 

rH 

t-4 

h- 

W U ^ 

w 


►-i x 


2 

L"* 















UJ 

UJ 



UJ UJ 

"5 

3 



o 

1- < 1- 1- 

K 


00 LU 

z 



-J 

on 

II 


OJ 


m 



in 


o 

QC 

o 

< • 

D O 

o 

*5 

« • 


2 


o 

o 

< z <t <1 

<t 1- 

Z JOH 

g? 

£D 

z 











a. 

O 

< 

< rH 

rH 


< < 


»-* 

o 

CL 

S XX 

2 

< 

LU CL 

X 

2 

X 

2 

< 

UJ 

II 

II 


ii 


II 

ti 


II 


rH 

-J 2 

2 ^ rH 

it 

z z 

«wr 


H* 

rH 

CD 

OC Cd OC oc 

cc 

X 

X X 

X 

-J 

-j 

UJ 

£X 

> 












-J w 

*— 



W w 

o 

Z 



o 

O rH O O 

o 

-J 

•HOC 

CO 

CD 

Z 1- o 

< 

< 


<r 


< 

< 


<r 

LL 

c 

<r u- 

LL C 

n 


LL U- 

<r 

o 

O 

a 

on 

LL U- LL 

LL 

<X 

o o 

VJ 

V. 




2l 

z 

2 


z 


2 

z 


z 

»— 4 

o 

o *-< 

»-4 Q 

a 


*-4 L-4 

<r 

o 

o 

Q 


hI CM If! O 


uoooooooo 


91 



00 140 J=1 # NN1 211700000 


APPENDIX C - Continued 


ooooocoooooo 

CJOOOOOOOOOOO 

OOOOOOOQOOOO 

CJOOOOOOOOOOO 

oooooooooooo 




o o o o o 
o o o o o 
o o o o o 
o o <_> o o 
o o o o o 

tn m ro m fn 


o 

o 

o 

o 

o 

in 

ro 


o 

o 

o 

o 

o 

o 

ro 


oooooooo 

OOOOOOOO 

OOOOOOOO 

oooooooo 

oooooooo 






















oo 
























U 

— » 








— . 















X 

</) 








CO 















» 

o 








ro 















1- 

X 






— • 


Z 















< 

» 






CO 


z 















* 

< 






CO 


•> 















t/> 

♦ 






z 


H 















o 

LO 






z 


II 















5“ 

u 






*» 


M 












— * 



•» 

X 






H 

















< 

► 






II 


0-^ 




•» 








»■ 



•* 

H- 






4—4 


CM 




4— 4 








M 



O 

< 






*> 


CM 















► 

•» 








Z 




< 








< 



rl 

o 






r—4 


z 















z 

•• 






Z 


•* 




II 








II 



z 

H 






z 


CM 










eg 






z 






*• 


Z 










cm 


•— > 



fH 

z 




z 


l“4 


z 


CM 


< 


o 




z 


< 



z 





a: 


II 


II 


CM 


"5 


o 



H 

z 





z 

H 

o 



3 


3 



z 

z 


•» 


eg 



z 

+ 


♦ 



« 

z 

H 



K 


♦ 


•» 

z 

z 


< 





z 



< 



</> 

z 




UJ 

CO 




* 

+ 


4-4 


o 



♦ 






z 

* 

o 



Od 

CO 

•*> 








h 



“5 

II 





<*• 

LU 

h- 




z 

♦ 



II 

M 


< 





ft 



< 



CL 

> 





z 

M 


►-4 

< 

< 


< 


O 



rH ^ 

< 


< 



h- 

o 

o 



o 

» 



w 






o 


z 

Z 1 

“5 






X 

o 



• m 

< 


< 




II 




z 

z 



II 



w* 





UJ 

rO 

<X 


<1 


— 






♦ 

» — 

— 

rg 




X 

X 




z 

z 




in 

H 

-*■ 





*-i m 

>C 

rg 

— - 



*— 4 

»-4 

fM 



• 

z 

w 


w 

• 

• 

z 

< 


• 

o 

11 

M • 

• z < 


o 

QC 

a: 

• 

O 


K 





a 

o 

z 

-> 

UJ 

a 

o 

•—I 

“5 O 

a z 

-MUO 

K- 

►— 

a o a* z 

• 

•> 


<► 

LU 

LU 

♦ 

•> 

3 uj 



UU 

UJ 

+ 

• 

3 ^ 

< 

<r 

LU 

H 

-J 

►—4 

H 

CM 

in 

CM 

• 

• 

4—4 

< 

z 

• 


O 

O “5 • 

a 


<r z 


X 

X 

• 


Z a: 






< 


►— < 

»-» < 

a >o 

<o < 

<L 


*-* 

NH o 



< 

a 


CL 


K 

h- 

H*- 

z 

z 

II 

w 

K 

z 

H- 

I— 4 

H II Z 

Z 

u 

w 

V- 

H* 


-j 

z 

K 

H- 

*•-4 

z 

Z 

Z 

Z 


w 


o 

Z 

w 



— ■ 

w 


o 

Z 


-J 

-i 

W 


Z 

w 


»— 4 

i-m 

►—< 

U- 

UL 

< 

< 

o u. 

O 

o 

O <r u. 

LL 

<r 

< 

c o 

< 

< 

LL 

O 

O U- 

cc 

CL 

CL 

ac 

*— * 

►-4 


< 

o 

►—4 

a o 

O ~3 >-M 

fc-H 

M 

< o o 

o 

o 

4— « 

o 

O 


a 

a 

Ql 

a 





o 



o 





o 


o 




O 













o 





o 


o 




o 










4— < 








^■4 


CM 




H 







92 


PRINT 5 214500000 
PRINT 2, ( ( AA( I» J) » J=NN23»NN33 ) » 1=1* NN33 I 214600000 
RETURN 214700000 
END 214800000 
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247 300000 

500 MOVE = 1 247400000 

CALL MATRIX (MOVE « l R0 ,IC0 t 0 ,0 1 MCS ,DUM ,MCS » 247500000 

RETURN 247600000 

END 247700000 



APPENDIX C — Continued 


ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
oo o o-H<vm«»-ir\vor~ t c 
r-h-ooaoocooooocoooco 


OOOoOOOOOO 

OOOOOOOOOO 

OOOOOOOOOO 

OOOOOOOOOO 

OOOOOOOOOO 

<X)QQ'Q'0'0'(?'0'Cr , 0' 

(N(\jcgrg<\JCV4<\jc\i«MC\i 


(\J 

» 

* 

3 


f\J 


m 

tt 

tt 

QC 

< 

CD 

V— 

* 

LU 


OO 

CL 

UJ 


* 

(VJ 

♦ 



Q 




00 






♦ 



•*• 









•> 





o 

♦ 



*■ 





WO 





JL 




o 

3 



3 





00 




'C 

U 




H 

X 



X 




X 

UJ 




H 

X 





1 



1 


<■*» 


U- 

z 




w 

Q 




a 

• 



• 


H 


o 

LL 




o 

Z 




K- 

f"-4 



H 


1 


a. 

ul 




» 

00 





w 





CD 


» 

►M 




3 

M 




o 

* 



X 


3 


X 

H 




X 

» 




o 

• 



oc 


wo 


«vr 

OO 




• 

K 





(M 



< 


z 


< 





UJ 

Z 




— » 




CD 


<w 


o 

O 

w 



» 

*— • 



— - 

rH 

w 



V— 


* 


z 

Z 

X 



2 

QC 



X 

• 

X 


X 

# 




3 


LL 



O 

a 



w 

o 

m 



• 


X 


LL 

o 

o 



X 

»— < 



00 

UJ 

* 


X 

Cvj 


* 



z 

U- 



CL 




00 

0 

* 


•— » 

X 


C 


UJ 

UJ 




• 

K 

cc 

CD 

UJ 

X 

CL 


X 

OsJ 




2 

GO 

z 



o 

CL 

< 

X 

z 

o 

< 

O H 

# 


II 




o 



z 

o 

CD 

O') 

X 

3C 

CD 

o 

CL 

* 

tu 



►— 



o 


<1 


f— 

Z 

X 

c 


(V 

UJ 

X 

3 



3 

UJ 

oo 3: 

z 




X 

H 

z 

* 


O 

* 

z 

r*^ 

z 

o 

»- 

z 


LJ 

m 

>0 

r- o' 


wo 

uj c 


Uj 

►-* 


CL 

a: 


UJ 

— 1 

X 

X 

X 

X 

X 

-J 



H- 

~J 


H* 

X 

3 

CD 

3 

X 

<1 

X 

-i 

-J 


—1 

— J 

V 

h 


«J 

1) 

z 

U- 

V- 

3 

— j 

*— 

UJ 

a 

CD 

CO 

CD 

CD 

<1 

X 


o 

<r 


a 

c 

LU 

00 


o 

a 

o v 

V. 


X 

o 

M 

o o o 

O 

o 

LL 

QC 


> 

UJ 


W H h (T 


o o 

a o 

rH C\| 


103 



APPENDIX C — Continued 


oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

oooooooooooo 

O'OrlNm^lfl lOMDtfO 
O'OOOOOOOOOO-H 
•TUMninirnninir.uurirMo 


OOOOOOOOOO _ 
OOOOOOOOOOO 
OOOOOOOOOOO 
OOOOOOOOOOO 
OOOOOOOOOOO 

H H H H H H H H H CM fSj 

in m in in in m if, m in in m 
cmcmcsicmcmcmcmcmcmcmcm 








oo 


















a 


















UJ 


















» 













z 





00 


















00 













<* 





•k 





o 








UJ 






X 




o 







— * 

oc 





-a 

o 




r— < 







X 

< 





r-* 









X 



LL 







o 




o 




* 



o 

f- 





a 

z 




h* 




LL 



u. 

*— • 





* 

oo 








O 



•» 

z 





x 

»*■« 




O 




X 



X 

3 





z 





o 




CC 


— • 

«— > 






♦ 

h- 








+ 


H 

00 

C* 



H 


UJ 

z 




— * 




CC 


1 

o 

UJ 





•» 

*—< 




rJ 



*■* 

< 


*—« 

z 

a 



X 


"X 

a: 



X 

• 



X 

CD 


w 

3 




u 


a 

a 




o 

a: 


•> 

K 


* 

LL 

00 



o 


X 

*«* 



00 UJ 

< 


X 

* 


l/> * 


00 



u 


(X 

• 


LL 

00 

• 

CO 


N»* 



• X 

UJ 

< 





♦ 

h- 

or. 

o 

UJ 

X 

h- 


X 

O 


*~4 * 

Z 

z 



z 


o 

Q- 

< X 

z 

o 

* 

O 


X 


II z 

*— « 




o 



O 

an 

CC 

X 

.A 

3S 

o a 

CC 

UJ 

H-. || 

h- 

CO 


3 

*— < 


< 


K 

X 

X 

ci a 

CM Uj 

* 

a 

«-* 


UJ 

z 

Z 

oo 

z 

X 

X 

X O 

h 

z 

X 


o 

* z 

o ►- 

a 

b- 



z 

o 

ro 

o 

h- in 


CO 

CC 

U 


C\J 

*— * 

o ~ 

cC 

< 


-J UJ 

z 

X 

* 



-J 

*— « 



— J 



cn x 

an 



< z 

z 

-J 

-J 

_i 

-J 

-J 

w 

(I 



1! 

Z 

UL 

3 

—I 

UJ 

UJ 

kw4 

o 

cr 

CD 

CC 

CD 

< 

LL 


o < 


o 

O o 

00 

< 

ac 

a: 

o 

u 

X 

X 

X 

X 

o 


z 

12 U 

Z 

u 

o u 


> 





m 

H 

h* 

H 










Ui 













o 


o 

o 















o 


o 

o 

















CM 



O 


104 



APPENDIX C — Continued 


OOOOOOOOOOOOOOOO 

oooooooooooooooo 

OOOOOOOOOOOOOOOO 

oooooooooooooooo 

oooooooooooooooo 

(MfNJfVIlMfNJPvICvlfMrrifriOfOOOOrri 

n->iriririinini^ir\iovriiAinir>inir. m 

cucmcmcmcucmcncmcmcmcmcmcmcmcmcm 


0 c 

Hi 



CO 




CM 






X 




* 






< 




* 






o 




*■* 










• 






o 




f\J 






z 




v. 










— 






•> 




44 






00 




o 






44 




♦ 










CM 






h- 




u 










44 






X 




1 






H 




X 










O 






1 


X 


o * 








o 


H 






U 


7 


44 






*«4 


a 


o u 






O' 


z 


K 1 





X 

K 


00 


CM 





4 

Ui 


44 


O U 





X 

X 


4 


O W 






X 


h- 


\ 





X 

>- 


Z 


— o 





K 

00 


44 


o * 

4^ 




Cl 



cm at 


• • 

X 




LU 

UJ 


o a 


o *■ 

44 




O 

eo 


4 44 


UJ 1 

oo 







44 4 


• **■ 

00 




UJ 

H 


O H 

Ct 

I -« 

UJ 




z 

00 


* a 

< 

u o 

z 




1-4 

3 


> a 

ou o 3 t 

X 

ac 



H* 

X 


'V. 44 

h* 

N U C\J 

X 

<i 



3 


z 

f*0 X 

^ ro z u z 

H* 

GO 

z 


O 

-J 

o 

t~4 sO 

r- _i »/) — a 


H- 

0£ 


QC 


X 

* X 

* * •- # 3 

—I 


3 


CD 

c 

X 

— J -J 

-J 

-1 — o H 

-J 

II 

h- 

a 

3 

U- 

o 

CO CD 

CD 

CC LL II UJ 

< 


UJ 

z 

00 

cc 

o 

V V 

N 

^ *-* X a: 

o 

X 

a 

UJ 


44 


44 r4 

h- 00 






< 





o 





o 


o 


105 
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end 256700000 
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The following is a listing of input data and a printout from a computer run used to determine the flutter- 
solution neighborhood: 
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CMEGA INITIAL CMIMT = 0. RADIANS/SEC 

CMEGA INCREMENT OELCM = A.0000000CE*00 RADIANS/SEC 

CMEGA FINAL CMFIN = A.OOOOOOOOE«-Ol RADIANS/SEC 



INITIAL V1NIT = 0. IN. /SEC 

INCREMENT CELV = 1 .00000000 E+OA IN. /SEC 
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The following printout is an optimization run with short printout option. The initial values of the flutter veloc- 
ity and flutter frequency are taken from the previous computer run. 
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CESIGN VARIABLES FINITE DIFFERENCE INCREMENTS ANALYTICAL DERIVATIVES * INDICATES ACTIVE 

DESIGN VARIABLES 
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APPENDIX C — Continued 
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APPENDIX C — Continued 
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APPENDIX D 


RESULTS FOR TRAPEZOIDAL FLAT PLATE 

In order to check the validity of the structural analysis, this appendix presents a 
comparison of theoretical and experimental natural frequencies of vibration of a canti- 
levered trapezoidal plate. The theoretical flutter speed and frequency for that plate are 
also presented but are not compared with experimental results. The dimensions and 
material properties of the aluminum plate are given in figure 12. The results are sum- 
marized in table I where theoretical vibration results are shown from SWIFT, which is 
the computer program based on the analysis presented in this report and described in 
appendix C, and from NASTRAN, a general-purpose finite- element structural-analysis 
computer program. The experimental results and the theoretical results obtained by 
using NASTRAN are taken from reference 26. 



1 

~T 

.250 in. (.635 cm) 

Figure 12.- Dimensions and material properties of aluminum 
plate used to obtain results given in table I. 

E = 1.00 x 10^ ksi (68.9 GK/m^) ; u = 0.30; Density, 

0.100 lbm/in? (2770 kg/m3). 

Achieving a true "fixed” boundary condition experimentally is very difficult. The 
support must be massive, and the specimen-attachment bolts must be numerous, heavy, 
and well-placed. The experimental results presented in table I are believed to approach 
a true fixed-boundary condition. 
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TABLE I.- EXPERIMENTAL AND THEORETICAL NATURAL FREQUENCIES 
AND THEORETICAL FLUTTER SPEED AND FREQUENCY OF 
CANTILEVERED TRAPEZOIDAL PLATE 


SHOWN IN FIGURE 11 


Experimental 
natural frequency, 
Hz 

Theoretical natural frequency, Hz, for — 

NASTRAN 
(109 grid points with 
327 degrees of freedom) 

SWIFT 

(10 spanwise stations with 
33 degrees of freedom) 

12.3 

12.1 

12.1 

38.6 

41.8 

42.9 

63.5 

63.7 

61.4 

95.0 

96.9 

107.0 

151.3 

147.7 

150.2 


Theoretical flutter speed . . 
Theoretical flutter frequency 


1180 m/s (3890 ft/s) 
35.7 Hz 









APPENDIX D - Continued 


A comparison between the experimental results and both sets of theoretical results 
indicates that the structural model used in SWIFT is adequate for simple wing structures. 

In order to provide some insight into the flutter calculations, curves showing the 
locus of points for which the real part of the flutter determinant is zero and curves 
showing the locus of points for which the imaginary part of the flutter determinant is zero 
are shown in figure 13. For this case, the two flutter solutions shown in figure 13 are 
widely separated. 

Following figure 13, a portion of the computer printout which contains the pertinent 
input data for the flutter calculations is presented. The input data used to obtain the 
natural frequencies differ only in that RHOA, the density of the air, was set equal to zero. 


Frequency, Hz 



Figure 13. - Locus of points for which the real part of the flutter 
determinant, R(Det), is equal to zero (solid lines) and for 
which the imaginary part, Im(Det), is zero (dashed lines). 
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DETERMINATION OF THICKNESS DISTRIBUTION FOR MINIMUM WEIGHT WING 
SATISFYING PRESCRIBED FLUTTER AND STRENGTH CRITERIA. 

REFERENCE - NACA TN-3640, STEIN- SANDER S> JUNE, 1956. 
OEXTER-STROUO RDK288 PAY, 1969. 
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OESIGN VARIABLES FINITE DIFFERENCE INCREMENTS ANALYTICAL DERIVATIVES * INDICATES ACTIVE 

DESIGN VARIJUU.ES 

Cl 1) = 2 . 50000000E-C 1 DC< II =0. ADI 1 1 = 1 .46385774E+02 



STATIONS 
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